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Chapter  1 
Introduction 


In  the  past  and  present,  solar  arrays  used  in  space  have  been  operating  at  low  voltage  levels, 
typically  at  28V.  Future  solar  arrays,  however,  are  being  designed  for  much  higher  voltages 
in  order  to  meet  high  power  demands  of  the  order  of  1  OkW  to  1  MW.  High  current  levels  could 
be  used  instead  to  achieve  these  increased  power  demands,  but  the  power  distribution  cables 
would  need  to  be  more  massive  and  the  resistive  losses  would  be  greater.  Consequently,  the 
current  is  maintained  at  a  low  value  while  the  voltage  is  increased  to  attain  the  necessary 
power  level. 

A  schematic  of  a  conventional  solar  cell  is  shown  in  Figure  1 .  The  coverglass  and  substrate 
shield  the  solar  cell  from  the  environment,  mainly  to  reduce  radiation  degradation.  These 
are  attached  to  the  cell  with  adhesives.  The  solar  cell  itself  is  a  semiconductor  of  two  parts, 
a  p-type  semiconductor  which  has  an  abundance  of  electrons  and  an  n-type  semiconductor 
which  has  an  abundance  of  electron  holes.  This  construction  allows  the  solar  cell  to  use 
the  photoelectric  effect  to  convert  solar  energy  into  electric  power.  A  photon  with  energy 
equal  to  or  greater  than  the  energy  gap  of  the  solar  cell  it  enters  will  free  an  electron.  This 
creates  an  electron-hole  pair.  If  the  pair  is  in  the  p-type  semiconductor,  the  electron  will  be 
accelerated  across  the  p-n  barrier  to  the  n-type  semiconductor  where  it  will  recombine.  The 
hole,  however,  will  be  repelled  by  the  barrier  because  of  the  excess  of  holes  in  the  n-type 
semiconductor.  Likewise,  if  the  electron-hole  pair  is  in  the  n-type  semiconductor,  the  hole 
will  be  accelerated  across  the  p-n  barrier  and  the  electron  repelled.  Consequently,  metal 
interconnectors  connect  the  n-type  semiconductor  of  one  cell  to  the  p-type  semiconductor 
of  the  adjacent  cell  to  utilize  the  current  created  by  the  electron  and  hole  movement.  Solar 
cells  are  connected  in  parallel  with  metal  interconnectors  to  obtain  desired  current  levels 
and  connected  in  series  to  obtain  desired  voltage  levels. 

The  solar  array,  along  with  other  surfaces  of  the  spacecraft  which  can  allow  the  passage  of 
current,  collects  current  from  the  ambient  plasma.  In  steady  state,  the  spacecraft  is  grounded 
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with  respect  to  the  plasma  by  the  zero  net  charging  condition 

dp  -! 

■^  +  Vj  =  0,  (1) 

which  is  derived  from  Ampere’s  Law  and  Gauss’s  Law.  To  obey  this  condition,  most  of  the 
solar  array  floats  negatively  with  respect  to  the  plasma.  This  is  because  the  random  thermal 
flux  of  the  lighter  electrons  to  the  spacecraft  is  greater  than  the  random  flux  of  the  heavier 
ions.  Therefore,  the  spacecraft  surfaces  must  be  at  a  negative  potential  in  order  to  maintain 
zero  net  current  collection. 

High-voltage  solar  arrays,  however,  have  been  observed  to  interact  with  the  plasma 
environment  of  low  earth  orbit  in  two  undesirable  manners.  For  positive  voltages,  the 
current  collection  can  be  anomalously  large,  possibly  leading  to  surface  damage^^  This 
phenomenon,  known  as  “snapover”,  occurs  when  the  dielectric  surface  potential  becomes 
positive,  attracting  electrons.  Above  a  certain  potential,  more  than  one  secondary  electron 
is  released  by  the  incident  electrons.  These  excessive  secondary  electrons  are  collected  by 
the  interconnector  and  seen  as  a  current  increase,  which  then  incurs  a  power  loss.  For  large 
negative  voltage  operation,  arc  discharges  can  occur".  Arcing  observed  in  experiments  has 
been  defined  as  a  sharp  current  pulse  much  larger  than  the  ambient  current  collection  which 
lasts  up  to  a  few  microseconds.  This  current  pulse  is  usually  accompanied  by  a  light  flash  at 
the  edge  of  the  solar  cell  coverglass.  Arc  discharges  can  cause  electromagnetic  interference 
and  solar  cell  damage^®,  so  there  is  a  need  to  study  mitigation  methods  and  to  be  able  to 
predict  arcing  rates  with  models. 

1.1  Background 

Arcing  has  been  studied  in  many  experiments  and  theoretical  arguments.  The  Plasma  Inter¬ 
actions  Experiments  (PIX I  and  II)  and  Solar  Array  Module  Plasma  Interactions  Experiments 
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(SAMPIE)  have  been  the  only  space  experiments  so  far,  though  several  more  space  experi¬ 
ments  are  planned  for  the  near  future.  Arcing  has  also  been  observed  in  many  ground  tests 
conducted  in  vacuum  plasma  chambers.  Two  different  theoretical  explanations  were  given 
by  Parks  et  al.^'  and  Hastings  et  al.'*’.  Cho  and  Hastings  ^  used  ideas  from  both  to  present 
a  more  complete  theory  of  the  arcing  sequence  of  events. 

Arcing  on  solar  cells  was  originally  observed  by  Heron  et  al.“  in  1971  during  a  high 
voltage  solar  array  test  in  a  plasma  chamber.  The  array  was  biased  to  -  16kV,  and  arcing  was 
observed  as  low  as  -6k V  in  a  plasma  density  of  10®m"^. 

In  1978,  the  first  Plasma  Interactive  Experiment  (PIX)®  confirmed  that  arcing  occurs  in 
space.  As  an  auxilliary  payload  on  Landsat  3’s  Delta  launch  vehicle,  PIX  operated  for  4 
hours  in  a  polar  orbit  around  920km.  A  solar  array  of  twenty-four  2cm  x 2cm  conventional 
silicon  cells  was  externally  biased  to  -lOOOV.  Arcing  discharges  began  at  -750V. 

In  1983,  PIX  IF  was  launched  also  as  an  auxilliary  payload  aboard  a  Delta  launch  vehicle 
into  a  near  circular  polar  orbit  of  approximately  900km  in  altitude.  The  500  2cm  x  2cm 
silicon  conventional  cells,  biased  to  -lOOOV,  experienced  arcing  as  low  as  -255V  and  at 
densities  as  low  as  lO^cm"^.  The  results  also  found  arcing  to  be  the  most  detrimental  effect 
of  negative  biasing. 

Ferguson  studied  the  PIX  II  ground  and  flight  results.  The  interconnectors  collected 
current  proportional  to  the  applied  voltage  bias.  The  arc  rate  R  was  determined  to  scale  as 

where  a  ~  5  for  the  ground  experiments  and  a  ~  3  for  the  flight  experiments,  ric  is  the 
ambient  plasma  number  density,  T,  is  the  ambient  ion  temperature,  and  rrii  is  the  ambient 
plasma  ion  mass.  The  dependence  of  the  arc  rate  on  these  parameters  indicates  that  the 
coverglass  surface  is  recharged  by  the  thermal  flux  of  ions. 

Ground  experiments  revealed  more  characteristics  of  the  arcing  phenomenon.  Experi¬ 
ments  by  Fujii  et  al.^  showed  that  dielectric  material  near  the  biased  conductor  in  the  plasma 
environment  is  essential  for  arcing  to  occur.  Fujii  et  al.  tested  material  plates  biased  to  high 
negative  voltages  in  a  plasma  environment.  The  plate  partially  covered  by  a  20{)//m  thick 
coverglass  experienced  arcing  at  -450V  while  the  uncovered  plate  did  not  arc,  except  at 
-lOOOV  when  the  arc  occurred  at  the  substrate.  Snyder^^  measured  the  electric  potential  on 
the  coverglass  and  found  that  it  decreased  significantly  when  an  arc  occurred.  This  indi¬ 
cates  that  the  negative  charge  created  during  arcing  discharged  the  positive  surface  charge 
accumulated  on  the  coverglass  surface.  Both  Snyder  and  Tyree^^  and  Inouye  and  Chaky'^’ 
observed  electron  emission  from  the  solar  array  that  could  not  be  explained  by  the  ambient 
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plasma.  Finally,  electromagnetic  waves  generated  from  the  arcing  current  were  measured 
by  Leung 

The  first  theoretical  model  was  proposed  by  Jongeward  et  al.'®  and  later  expanded  by 
Parks  et  al.^^ .  Jongeward  et  al.  attributed  Snyder’s^^  experimental  observation  of  the  decrease 
in  coverglass  potential  prior  to  arcing  to  enhanced  electron  emission  from  the  interconnec¬ 
tor,  which  corresponds  to  the  electron  emission  observed  in  References  23  and  15.  They 
suggested  the  emission  is  due  to  a  thin  layer  of  ions  deposited  on  the  interconnector,  causing 
the  electric  field  to  be  significantly  increased.  The  time  for  positive  charge  build  up  is  then 
dependent  on  the  ambient  density  n^,  the  interconnector  size,  and  the  bias  voltage.  The 
arc  discharge  is  proposed  to  occur  by  a  positive  feedback  mechanism  from  electron  heating 
which  leads  to  a  space  charge  limited  discharge.  At  low  ion  densities,  other  surface  neutral¬ 
izing  effects  are  said  to  dominate,  thus  inhibiting  the  positive  charge  build  up.  Jongeward  et 
al.  also  modeled  the  arc  discharge  decay  time  by  assuming  space  charge  limited  conditions 
and  showed  that  the  peak  current  magnitude  agrees  well  with  this  assumption. 

Parks  et  al.^^  concentrated  on  further  detailing  the  theory  proposed  by  Jongeward  et 
al.’®  on  the  prebreakdown  electron  emission  current.  They  accepted  Jongeward’s  theory 
of  positive  charge  build-up  in  a  thin  insulating  layer  on  the  interconnector  and  of  arcing 
orginated  from  interconnector  electron  emission  instead  of  from  the  ambient  plasma.  Parks 
et  al.  proposed  the  addition  of  the  phenomena  presented  by  Latham”'  ’®,  namely  that  non- 
metallic  emission  processes  are  significantly  responsible  for  electron  emission  by  nominally 
metallic  surfaces.  Therefore,  Parks  et  al.  claimed  that  the  arc  rate  must  be  proportional  to  the 
electron  emission  current  density  and  the  bias  voltage.  They  further  suggested  that  electron 
emission  is  controlled  by  the  vacuum  electric  field  at  the  surface  of  the  insulator.  Given 
these  assumptions,  they  determined  that  the  rate  of  field  build-up  in  the  insulator  is 

—  Exns—vac)  —  ~  1),  (3) 

where  tins  is  the  dielectric  constant  of  the  insulator  layer,  E^s  is  the  electric  field  inside  the 
insulator,  Ems-vac-  is  the  electric  field  at  the  insulator- vacuum  interface,  j,  is  the  ion  current 
density,  jpN  is  the  Fowler-Nordheim  emission  current  at  the  metal-insulator  interface,  a  is 
the  rate  of  ionization  per  unit  distance  inside  the  layer,  d  is  the  thickness  of  the  insulator 
layer,  and  Pis  the  probability  that  electrons  are  emitted  from  the  insulator- vacuum  interface. 
The  emission  current  from  the  metal-insulator  interface  is  given  by 

jFN  =  AElx^e  ,  (4) 

where  A  and  B  are  the  Fowler-Nordheim  emission  coefficients,  given  in  Eq.  (8)  and  (9).  This 
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Figure  2:  Model  of  the  conventional  solar  array  used  for  numerical  simulations 


expression  for  the  electric  field  accounted  for  experimental  observations  of  the  characteristics 
of  the  voltage  threshold,  the  prebreakdown  electron  emission  current,  and  the  arcing  rate. 

Hastings  et  al.'°  did  not  try  to  explain  the  prebreakdown  electron  emission  current  but 
instead  proposed  a  model  for  the  gas  breakdown  seen  as  the  arc  discharge.  They  suggested 
that  neutral  gas  is  desorbed  from  the  sides  of  the  coverglass  by  electron  bombardment,  a 
phenomena  known  as  electron  stimulated  desorption  (ESD).  The  bombarding  electrons  are 
emitted  from  the  interconnector,  as  determined  from  Snyder’s  experiments^^,  and  from  the 
coverglass  as  secondary  electrons  which  return  to  the  side  surface.  The  desorbed  neutrals 
then  accumulate  in  the  gap  between  the  coverglasses  over  the  interconnector,  forming  a 
potentially  high-pressure  gas  layer  which  can  break  down  from  the  electron  emission  current 
flowing  through  it.  This  was  in  contrast  to  the  previous  theory  which  suggested  that  the  arc 
occurs  in  an  insulator  on  the  surface  of  the  interconnector. 

Recent  work  by  Cho  and  Hastings^  combined  some  of  the  ideas  from  these  two  theories 
and  studied  the  charging  of  the  region  near  the  plasma,  dielectric,  and  conductor  triple 
junction.  The  model  that  they  studied  is  shown  in  Figure  2.  The  dielectric  consists  of  both 
the  coverglass  and  the  adhesive  bonding  the  coverglass  to  the  solar  cell.  The  conductor  is 
the  interconnector,  which  is  usually  placed  between  the  cell  and  substrate  on  one  end  and 
between  the  cell  and  cover  adhesive  in  the  adjoining  cell.  The  solar  cell  itself  was  neglected 
since  the  potential  drop  across  it  is  at  most  a  few  volts  while  the  potential  drop  across  the 
coverglass  and  adhesive  is  hundreds  or  even  thousands  of  volts  for  high  voltage  operation. 

Cho  and  Hastings  developed  a  numerical  simulation  of  the  arc  initiation  processes.  They 
studied  charging  of  the  dielectric  surfaces  by  three  sources:  ambient  ions,  ion-induced 
secondary  electrons,  and  enhanced  field  electron  emission.  From  numerical  results,  they 
determined  the  following  arc  sequence,  illustrated  in  Figure  3: 

(1)  ambient  ions  charge  the  dielecU'ic  front  surface,  but  leave  the  side  surface  effectively 
uncharged; 

(2)  ambient  ions  induce  secondary  electrons  from  the  conductor  which  charge  the  side 
surface  to  a  steady  state  unless  enhanced  field  electron  emission  (EFEE)  becomes 
significant; 
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Figure  3:  Arcing  sequence  of  a  high  voltage  solar  array 

(3)  EFEE  will  charge  the  side  surface  if  there  is  an  electron  emission  site  close  to  the  triple 
junction  with  a  high  field  enhancement  factor,  /?;  and 

(4)  EFEE  can  result  in  collisional  ionization  of  neutrals  desorbed  from  the  coverglass, 
which  is  what  is  observed  as  the  arc  discharge. 

They  also  found  that  the  electric  field  at  the  triple  junction  is  not  bounded  during  EFEE 
charging. 

Cho  and  Hastings  used  the  numerical  results  to  develop  analytical  formulas  describing 
the  arcing  rate^’^  They  suggested  that  the  time  between  arcs  is  the  minimum  of  the  sum 
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of  the  ambient  charging  time  Ticm  and  the  enhanced  field  electron  emission  charging  time 
Te/ee,  SO  that  the  arc  rate  R  is  given  by: 

R  =  min{Tian  +  Te/eer^  (5) 


For  the  ion  charging,  Cho  and  Hastings  showed  that  ambient  ions  mainly  charge  the  front 
surfaces  of  the  coverglass,  not  the  side  surfaces.  They  expressed  this  time  as 


‘T'ion 


AQ 

CTlf  VicmAccii 


(6) 


where  AQ  is  the  charge  lost  by  tHe  coverglass  due  to  an  arc  discharge,  which  must  be 
recovered;  ericVion  is  the  ambient  ion  flux  to  the  front  surfaces,  with  Vion  as  the  mean  speed 
of  ions  entering  the  sheath  surrounding  the  solar  array;  and  Aceii  is  the  frontal  area  of  the 
coverglass.  Assuming  a  constant  secondary  electron  yield  and  constant  voltage  bias,  they 
derived  the  following  analytical  expression  for  the  EFEE  charging  time,  Tejee.'- 


c dieted}  (  Bd  \ 

"  (7„  -  \l3Vr,iJ 


where  A  and  B  are  the  Fowler-Nordheim  coefficients  given  by 


A  = 


1.54  X 

(pw 


(8) 


B  =  6.53  X  10^ 


(9) 


^  is  a  factor  to  account  for  the  difference  in  dielectric  constants  and  is  given  by 


—  ^<^2  {dj 

di  (-di  di 


(10) 


.^0  is  {  evaluated  with  cf,  =  d,  d,  is  the  distance  from  the  triple  junction  of  the  first  impact 
by  an  electron  emitted  from  the  conductor,  Cdteic  is  the  capacitance  of  the  dielectric  at  this 
impace  site,  7ee  is  the  secondary  electron  yield  of  the  dielectric,  Sreai  is  the  real  area  of  the 
emission  site,  Sfn  is  the  Fowler-Nordheim  “effective”  area  of  the  emission  site,  77  is  a  factor 
to  account  for  the  difference  in  the  electric  field  at  the  emission  site  from  the  triple  junction, 
d]  is  the  thickness  of  the  coverglass,  <^2  is  the  thickness  of  the  adhesive,  d  =  di  +  d2,V  is  the 
voltage  at  which  the  arc  occurs,  and  P  is  the  field  enhancement  factor.  From  comparison 
with  experiments,  Hastings  et  al.®  suggested  that  P  must  be  greater  than  a  few  hundred,  so 
they  assumed  that  the  field  enhancement  is  due  to  a  thin  dielectric  layer  on  the  conductor 
surface  rather  than  microprotrusions.  They  later  updated  their  views  as  explained  in  Section 
2.1.1. 
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From  experimental  observations,  Hastings  et  al.®  suggested  other  characteristics  of  the 
arcing  processes,  such  as  the  discharge  wave  hypothesis  and  the  occurrence  of  one  arc  at  a 
time  within  a  certain  area.  The  discharge  wave  hypothesis  suggested  that  at  arc  initiation 
emitted  electrons  form  a  plasma  cloud  over  the  solar  array.  Some  of  the  electrons,  attracted 
by  the  positive  surface  potential,  strike  the  coverglasses  in  the  local  area  until  they  are 
discharged.  Experimental  results  also  showed  that  the  arc  current  is  more  likely  to  be  carried 
by  electrons,  consistent  with  the  hypothesis  that  arcing  is  initiated  by  electron  emission  from 
the  interconnector.  In  addition,  as  the  temperature  increased  fewer  neutral  gas  molecules 
were  desorbed  from  the  dielectrics  and  the  arc  rate  was  seen  to  decrease,  consistent  with  the 
hypothesis  that  ionization  of  the  neutral  gases  also  plays  a  role  in  arc  initiation. 

1.2  Overview  of  This  Research 

Power  requirements  for  space  systems  are  increasing  significantly.  As  the  most  reliable 
power  source,  high  voltage  solar  arrays  will  be  needed  to  meet  these  requirements.  Since 
arcing  degrades  the  array  performance  and  causes  electromagnetic  interference  which  affects 
nearby  instruments,  it  is  imperative  to  study  arcing.  Recent  studies  by  Cho  and  Hastings^’^ 
determined  an  arcing  sequence  of  events  and  an  arcing  rate  based  on  numerical  and  the¬ 
oretical  work  which  has  been  shown  to  agree  well  with  experimental  results.  With  these 
models  it  is  possible  to  determine  methods  of  arc  rate  mitigation  and  to  predict  arc  rates 
for  experiments.  This  research  can  then  be  used  in  the  design  of  new  solar  cells  and  in  the 
design  of  high  voltage  solar  arrays. 

The  focus  of  this  research  is  twofold:  to  identify  and  study  mitigating  effects  on  arc  rates 
and  to  present  arc  rate  predictions  for  the  Air  Force’s  Photovoltaic  Array  Space  Power  Plus 
Diagnostics  (PASP  Plus)  experiment.  In  Chapter  2,  the  numerical  and  analytical  nodels 
developed  by  Cho  and  Hastings  are  reviewed,  and  the  numerical  model  modified  for  the 
wrap-through-contact  solar  cell  geometry  is  presented.  Based  on  the  analytical  model  for 
conventional  solar  cells,  arc  rate  reduction  methods  are  studied  using  the  corresponding 
numerical  model  in  Chapter  3.  In  Chapter  4,  arcing  rates  are  predicted  for  the  high-voltage 
biased  arrays  of  the  PASP  Plus  experiment.  Finally,  conclusions  are  summarized  in  Chapter 
5. 
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Chapter  2 

Numerical  and  Analytical  Models 


2.1  Conventional  Solar  Cells 

A  schematic  of  a  conventional  solar  cell  is  shown  in  Figure  1 .  In  high  voltage  operation,  the 
voltage  differential  over  the  coverglass  and  adhesive  can  be  hundreds  or  even  thousands  of 
volts  while  the  voltage  differential  over  the  cell  itself  is  at  most  a  few  volts.  For  modeling 
purposes,  the  cell  semiconductor  can  therefore  be  neglected,  as  shown  in  Figure  2.  In  this 
model,  the  interconnector  is  a  conductor  and  the  coverglass  and  adhesive  are  dielectrics. 
The  numerical  and  analytical  models  used  for  conventional  cells  were  developed  by  Cho 
and  Hastings^’®,  as  briefly  described  in  Section  1.2. 

2.1.1  Numerical  Model 

The  numerical  model  incorporates  all  relevant  physical  characteristics  and  processes  for 
solar  cell  charging  from  the  ambient  plasma,  electron  emission  from  the  interconnecior, 
and  secondary  electron  emission  from  the  dielecuics.  A  representation  of  this  system  is 
shown  in  Figure  4.  The  coverglass  and  adhesive  surface  charge  densities  are  affected  by 
the  ion  ram  current  density  jid,  the  electron  emission  current  density  from  the  conductor 
jec,  and  the  secondary  electron  current  density  from  the  surface  jee.  After  arc  initiation,  the 
current  densities  from  the  ionization  of  neutral  gases  may  also  be  significant.  These  are  not 
considered,  however,  as  only  the  time  to  arc  initiation  is  the  focus  of  this  research.  The  rate 
of  change  of  the  dielectric  surface  charge  density  can  then  be  expressed  as 

=  jtdixA)  -  J  P{x,y,t)jec{y,t)dy  -  J  P{x,x\t)j^c{x')dx'  +  j.A'x.t)-  (11) 

where  P(x,y,t)  is  the  probability  that  an  electron  emitted  from  position  y  on  the  conductor 
hit  the  dielectric  at  position  x  at  time  t,  and  P{x,x',t)  is  the  probability  that  an  electron 
emitted  at  x'  hits  the  dielectric  at  x. 
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Figure  4:  Model  system  of  the  high  voltage  solar  array  and  plasma  interactions 

The  numerical  model  consists  of  three  schemes.  The  first  scheme  uses  the  capacitance 
matrix  method  to  obtain  a  preliminary  electric  potential  distribution  along  the  dielectric 
surfaces.  The  second  scheme  involves  a  particle-in-cell  (PIC)  method  which  is  used  for 
ambient  ion  charging.  Once  a  steady  state  is  obtained  from  the  ion  charging,  a  space- 
charge-free  orbit  integration  scheme  calculates  the  electron  charging  by  enhanced  field 
electron  emission  (EFEE). 

All  schemes  use  the  same  computational  domain  and  grid.  The  phase  space  of  the  domain 
consists  of  two  position  coordinates  and  three  velocity  coordinates.  As  shown  in  Figure  5, 
the  domain  includes  two  halves  of  solar  cells  with  the  interconnector  forming  the  lower 
boundary  of  the  gap  between  the  cells.  The  boundary  condition  far  from  the  cells  at  x  =  0  is 
0  =  0,  simulating  the  far  field.  In  these  simulations,  any  electrons  leaving  the  domain  at  x  = 
0  will  also  leave  the  sheath.  The  boundaries  are  thus  Dirichlet  in  the  x  direction  and  periodic 
in  the  y  direction  to  simulate  a  solar  array.  The  grid  is  clustered  along  the  dielectric  sides 
and  near  the  interconnector  for  better  resolution  of  the  large  electric  potential  gradients  in 
these  areas. 

The  capacitance  matrix  method  is  used  to  obtain  an  inititial  condition  for  the  PIC  code, 
thus  reducing  the  simulation  time.  In  employing  this  method,  which  is  given  in  Reference 
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14,  a  unit  charge  is  ascribed  to  one  cell  on  the  dielectric  surface  while  all  other  cells  have 
zero  charge.  The  Poisson  equation  is  then  solved  to  determine  the  electric  potential  in  every 
cell  on  the  dielectric  surface  due  to  this  unit  charge.  This  process  is  repeated  for  each  grid 
cell  along  the  surface  of  the  dielectrics.  Afterwards,  the  array  containing  all  of  the  potential 
values  calculated  is  inverted  to  determine  the  capacitance  value  for  each  grid  cell.  This 
matrix  is  stored  for  use  as  the  initial  conditions  of  the  PIC  code,  so  that  the  simulation 
can  be  started  from  any  charging  state  described  by  only  the  surface  charge  or  the  surface 
potential. 

With  the  capacitance  matrix  calculated  for  unit  charges,  the  PIC  code  calculates  the 
space  potential  based  on  a  pre-determined  surface  potential.  The  initial  conditions  for  the 
dielectric  potential  are  0  =  0  on  the  front  surface  and  a  linear  distribution  of  0  on  the  side 
surface  with  the  conductor  voltage  at  one  end  and  the  front  surface  zero  voltage  at  the  other 
end.  The  ram  velocity  is  oriented  90"  to  the  dielectric  front  surfaces  and  conductor.  To 
save  computational  time,  an  artificial  ion  mass  is  used  such  that  m^/me  =  100.  Ions  and 
electrons  are  initially  inserted  uniformly  throughout  the  domain  according  to  the  ambient 
density.  After  the  space  potential  is  calculated  using  the  Poisson  equation,  the  ions  and 
electrons  are  moved  according  to  the  new  potential.  A  new  space  charge  density  for  each 
grid  point  can  then  be  calculated  based  on  the  new  ion  and  electron  positions.  This  loop  is 
then  repeated  with  the  potential  being  re-calculated  based  on  the  new  charge  density.  The 
PIC  code  is  run  for  a  time  equivalent  to  the  inverse  of  the  ion  plasma  frequency  to  adjust 
the  space  charge  completely  with  the  surface  potential. 

The  results  from  the  PIC  scheme  are  the  initial  conditions  for  the  dielectric  charging 
scheme.  A  typical  contour  plot  of  the  initial  electric  potential  is  shown  in  Figure  6,  and  the 
corresponding  surface  charge  density  is  shown  in  Figure  7. 

No  electron  emissions  from  the  conductor  or  dielecu-ic  are  taken  into  account  in  the  PIC 
code  since  they  are  negligible.  The  electron  emission  which  leads  to  arc  initiation  was 
determined  to  be  enhanced  field  electron  emission  (EFEE)  by  Cho  and  Hastings^  They 
described  this  current  density  from  a  finite  emission  site  on  the  conductor  surface  as 

which  is  the  Fowler-Nordheim  expression  for  field  emission  due  to  a  thin  dielectric  layer  with 
the  added  factor  Spp^  / Sreai  to  account  for  the  negative  space  charge  effect  near  the  emission 
site.  The  electric  field  E  in  this  expression  is  the  electric  field  at  the  dielectric-vacuum 
interface.  A  and  5  are  the  Fowler-Nordheim  emission  coefficients  given  by  Eq.  (8)  and  (9). 
The  field  enhancement  factor  P  is  assigned  to  the  emission  site  to  represent  an  enhancement 
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Figure  6:  Typical  electric  potential  contour  plot  for  ambient  ion  charging  of  conventional  cells 
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Figure  7:  Typical  surface  charge  density  along  side  dielecU'ics  after  ambient  ion  charging  of  conventional  cells 
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Figure  8:  Electric  field  lines  over  a  whisker  on  conductor  surface 


due  to  manufacturing  defects  or  impurities.  As  shown  in  Figure  8,  the  protrusion  causes  a 
higher  electric  field  gradient  which  enhances  the  electric  field  at  its  tip.  From  electrostatic 
theory  for  a  whisker,  0  is  the  factor  of  the  enhanced  electric  field  at  the  tip  of  the  whisker 
defect  relative  to  the  average  electric  field  in  the  vicinity  is  equivalent  to  the  ratio  of  the 
height  of  the  protrusion  to  its  radius  of  curvature.  This  factor  can  be  of  the  order  of  1  to 
10^,  with  typical  values  of  interest  in  the  hundreds.  The  secondary  electron  current  density 
at  each  point  x  is  given  by 

Jee{x,t)=  J  %e{X,y)PiX,y,t)jec{y,t)dy-¥  J  Jee{X,x')P{X,x/,t)je.e{x')dx'.  (13) 

In  the  orbit  integration  scheme,  the  first  term  in  Eq.  (1 1)  is  neglected  since  it  was  shown 
to  be  insignificant  during  electron  charging^.  Using  Eq.  (13),  the  surface  charge  density 
equation  can  be  rewritten  as 

=  f  ilee{X,y)-\)P{X,y,t)j,,c[y0)dy-\-  J{'ree{X,x')-l)P{X,x\t)jee.{x')dx'.  (14) 

The  orbit  integration  scheme  then  consists  of 

(1)  obtaining  the  surface  potential  by  using  the  capacitance  matrix  method; 

(2)  solving  Laplace’s  equation  to  obtain  the  space  potential; 

(3)  integrating  test  electron  orbits  from  the  conductor  to  calculate  7ef  and  the  impact 
probabilities  P  for  a  given  electron  current  density  from  the  conductor; 

(4)  solving  Eq.  (1 1)  and  (13)  for  the  secondary  electron  current  density  j^r  and  the  rate  of 
change  of  the  surface  charge  density; 

(5)  renewing  the  surface  charge  density; 

(6)  obtaining  the  new  potential  for  the  renewed  surface  charge  density; 

(7)  calculating  the  timestep; 

(8)  determining  if  the  space  charge  current  density  is  too  high  or  the  timestep  is  too  small, 
either  of  which  will  halt  the  program;  and 
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(9)  calculating  electron  trajectories. 

Steps  (4)  through  (9)  are  repeated  until  the  specified  number  of  timesteps  are  completed  or 
the  progam  is  stopped  in  step  (8).  If  the  space  charge  current  density  is  too  high,  the  space 
charge  effects  of  the  emission  current  can  no  longer  be  neglected  so  the  PIC  code  must  be 
run  if  further  calculations  are  needed.  If  the  timestep  is  too  small,  the  electric  field  is  most 
likely  running  away. 

The  timestep  for  EFEE  charging  is  calculated  based  on  the  rate  of  change  of  the  dielectric 
surface  charge  density  at  the  first  impact  point  x  =  di.  This  can  be  expressed  by  neglecting 
the  second  term  in  Eq.  (11),  reducing  the  equation  to 

(7ee  -  l)yec(y,^)f^y-  (15) 

The  integral  P{x,y,  t)dx  is  approximately  unity  since  the  point  x^  d^is  the  first  impact 
point  by  emitted  electrons.  The  equation  then  simplifies  to 

=  (7ef  -  (16) 


or 

A(  = - — - 7=—,  (17) 

(Tee  1  )jec(y)  f  )(  v3'/dt) 

where  S  is  the  area  of  the  emission  site,  as  shown  in  Figure  9.  The  potential  difference 
between  the  triple  junction  (x  =  0)  and  the  impact  point  {x  =  di)  can  be  expressed  as 

=  (18) 
^diele 

where  CdteU  is  the  capacitance  of  the  dielectric  surface  and  a  is  the  surface  charge  density. 
The  electric  field,  then,  is  approximately 


E  = 


a 


^diclcdi 


The  timestep  can,  therefore,  be  determined  by  solving 


At  =  0.02 - - 

(Tee 


A  E  E  diclc  dt 
i)jec{yJ){E^/di)' 


(19) 


(20) 


where  the  empirical  factor  0.02  is  used  so  that  the  timesteps  will  be  shorter  than  the  actual 
timescale  of  arc  initiation. 


2.1.2  Analytical  Model 

The  analytical  model,  which  is  used  to  calculate  the  arc  rates,  is  drawn  from  the  theory  of 
Cho  and  Hastings^  ®,  discussed  in  Section  1.2.  The  arc  rate  is  determined  by  calculating  the 
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conductor 


Figure  9:  Geometry  for  EFEE  charging 


lime  between  arcs,  Tare,  given  by 

Tare  =  Tn^n{T^on  +  Te  fee),  (21) 

where  is  the  ambient  ion  charging  time  given  by  Eq.  (6)  and  Tefee  is  the  enhanced  field 
electron  emission  charging  time  given  by  Eq.  (7).  The  analytical  expression  for  Tefee  is 
determined  by  starting  with  Eq.  (16).  A  schematic  of  the  geometry  considered  is  shown  in 
Figure  9. 

The  electric  field  at  the  triple  junction  can  be  expressed  as 

<7 


^  Cameled,' 


(22) 


where  ^  is  given  by  Eq.  (10)  if  the  first  impact  site  of  the  electrons  emitted  from  the 
interconnector  is  on  the  coverglass  side  surface.  If  it  is  on  the  adhesive  side  surface,  E,  is 
unity.  The  electric  field  at  the  emission  site,  Ee,  can  be  very  different  from  the  electric  field 
at  the  triple  junction.  To  account  for  this,  the  factor  rj  is  introduced  so  that 

a 


Ee  =  i]Etj  =  . 

Substituting  Eq.  (12)  and  (23)  into  Eq.  (16)  results  in 

^ )  \/^real 


dEr 


dt  Celeeled^ 

This  can  be  integrated,  assuming  the  secondary  electron  yield  is  constant,  to  obtain 

Ee. 


he. 


-A'p^E^exp 


B 

'0Ee 


EeXt) 


1  +  (i  -  exp(-^^)(|)C'f  j  ’ 


(23) 


(24) 


(25) 
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Figure  10;  Typical  electric  field  nui-away  versus  time 


where  C  is  the  constant  given  by 


c  =  ni 


{'Ifee  ^  )  y/^real 
Cdieled^ 


A'0^ 


(26) 


and  is  the  initial  electric  field  at  the  electron  emission  site  on  the  interconnector.  From 
the  numerical  simulations  it  is  known  that  the  electric  field  usually  exhibits  the  behavior 
shown  in  Figure  10.  The  field  run-away  to  infinity  corresponds  to  the  denominator  of 
Eq.  (25)  reaching  zero.  This  run-away  time  also  corresponds  to  the  time  Xe/ee,  so  Xe/ee  can 
now  be  determined: 


'^e  fee  — 


l-exp(-^)  /fix 

‘exp(-^)  V/?y 

0  (  B  \\ 

c 

_ Cdreled'i _ ^  ^ 

(lee  -  ^)VSrealV^A§^^B0  ^KpE^.J 


(27) 

(28) 
(29) 


This  equation  is  the  same  as  Exj.  (7)  with  Ec„  expressed  as  the  potential  difference  between 
the  coverglass  front  surface  and  the  triple  junction; 

E,„  =  (30) 


Finding  the  minimum  of  the  sum  of  the  ion  and  EFEE  charging  times  accounts  for  the 
fact  that  EFEE  charging  can  initiate  whenever  the  surface  has  a  strong  enough  electric  field, 
not  just  when  the  front  surface  current  returns  to  zero.  To  find  this  minimum  charging 
time,  electron  emission  sites  must  be  considered  along  the  entire  conductor  as  opposed  to 
the  numerical  model  which  sets  one  emission  site  usually  next  to  the  triple  junction.  For 
each  emission  site  a  voltage  14  can  be  calculated  at  which  the  arc  occurred  by  solving  the 
differential  equation 
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or 


d 

dV, 


(K  -  (v;.-  -  ^ 


c 


frorx 


:))C, 


ront 


€-Tl^  Viofi 


+ 


cell 


Cdieledj^  f  Bd 

(7ee  -  ^)\/^real'niA§^B0  V^V^T/^O 

‘^real 


0. 

(32) 


where  14rc  is  the  voltage  of  the  last  arc  discharge  and  C front  is  the  capacitance  of  the 
coverglass  front  surface. 

In  order  to  solve  Eq.  (32),  a  number  of  properties  must  be  known  or  determined.  First, 
the  following  cell  properties  must  be  known:  the  thickness  of  the  coverglass  and  cover 
adhesive  (^1,^2),  the  dielectric  constants  of  the  coverglass  and  adhesive  {td,  ,td^),  the  energy 
of  incident  electrons  for  maximum  secondary  electron  yield  for  the  coverglass  and  adhesive 
(Bmaxi  ^^rnax2  ),  the  maximum  secondary  electron  yield  at  normal  incidence  for  the  adhesive 


and  coverglass  h’max,,lmax2),  the  interconnector  work  function  (0,^),  and  the  solar  cell 
frontal  area  {Aceii).  Then,  the  following  factors  can  be  determined:  A  according  to  Eq.  (8); 
B  from  Eq.  (9);  ^  and  from  Eq.  (10);  d  =  di  +  dj;  C/ront  which  is  approximated  as 


C 


front  — 


{Acclltd0/di  {Acclltd^) / d2 

and  7ef,  which  is  given  by  Reference  8 


(33) 


'fee.  —  ')max ' 


E. 


exp  [2  -  2\ 


E. 


exp[2(l  -  COS0,)]. 


(34) 


Here,  Ei  is  the  incident  energy  of  the  emitted  electrons  impacting  the  dielectrics  given  by 


E,  =  e(t)d  =  ^^  = 

^  C  d 


(35) 


and  6*1  is  the  incident  angle  of  those  electrons  at  the  first  impact  site  given  by 


d,  =  arctan 


(36) 


where  y  is  the  distance  of  the  emission  site  from  the  triple  junction.  The  mission  parameters 
determine  the  ion  velocity  v^on  and  the  range  of  the  ambient  density  n,.  If  the  array  is 
orientated  at  90°  to  the  ram  velocity,  v,oa  is  the  orbital  velocity.  Otherwise,  iw  is  a  sum  of 
the  orbital  velocity  and  the  mean  thermal  speed  of  ions  Ci/A,  where 


c,  = 


(37) 


Consequently,  the  ion  mass  and  electron  temperature  must  also  be  known.  For  each  arc 
calculation,  ric  is  chosen  randomly  from  a  uniform  distribution  in  log]o  rif .  Other  properties 
only  known  within  a  range  include  areas  Sefj  and  Sreai  and  enhancement  factor  0.  Areas 
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Seff  and  Sreai  ^rc  randomly  chosen  from  uniform  distributions  in  log,o  S^ff  and  logjo  Sreai, 
respectively,  between  given  minimum  and  maximum  values.  The  enhancement  factor  /?  is 
randomly  selected  from  the  distribution  f{P)  =  /oexp(-/?//?o),  where  /<,  is  determined  from 
the  normalization:  /  f{P)d(5  -  1.  Finally,  the  three  parameters  left  to  be  determined  are 
Cdieie,  dt,  and  rj,  all  of  which  are  functions  only  of  the  emission  site  distance  y  from  the  triple 
junction.  To  determine  Cdieie,  the  capacitance  matrix  scheme  used  with  unit  surface  charge 
values  must  be  run.  The  relevant  values  are  the  diagonal  elements.  Those  that  correspond 
to  the  lower  side  dielectric  are  non-dimensionalized  by  the  normal  capacitance 


C 

'^norm 


(38) 


and  inverted.  The  corresponding  distances  from  the  triple  junction  are  non-dimensionalized 
by  the  thickness  of  the  two  dielectrics,  d.  These  values  are  plotted  and  fit  to  a  five  order 
polynomial: 

n.=5  /  \  ^ 

=  .  (39) 

n=0  ^  ^ 

To  determine  di  and  r?,  results  from  the  orbit  integration  scheme  of  the  numerical  model  are 
used  to  obtain  functional  forms.  These  are 


d 


71—4 


(40) 


and 


n=4 


V 


(41) 


where  y  =  y/(dgap/2). 

The  voltage  Ve  is  determined  to  be  in  the  range  of  Vj,  the  voltage  differential  between 
the  front  surface  and  conductor  just  after  the  arc,  and  Vinas-  If  Te/er  dominates  to  the  point 
where  Tion  is  insignificant,  14  =  K-  Likewise,  if  r,o7i  dominates,  14  =  Vbias-  Otherwise,  the 
arcing  time  is  affected  by  both  Tf-jee  and  so  14  is  determined  by  the  Newton-Raphson 
method.  After  Tare  is  calculated  for  every  emission  site,  typically  numbering  1000,  the 
smallest  Tare  is  compared  with  the  experiment  time,  Texp-  If  Teip  is  greater  than  Tare,  another 
Tare  is  Calculated  until  the  sum  of  the  arcing  times  is  greater  than  Texp-  The  arc  rate  is  then 
the  number  of  arcs  counted  less  one  divided  by  the  experiment  time. 

For  a  given  solar  array,  the  surface  is  divided  into  sections  of  area  equivalent  to  the  area 
covered  by  the  arc  discharge  wave.  Based  on  experimental  measurements  in  Reference  9, 
this  area  is  chosen  to  be  0.0 12m^.  All  arcs  in  a  section  are  assumed  to  be  correlated,  but 
arcs  are  assumed  to  be  uncorrelated  between  different  sections.  The  arc  rate  is  calculated 
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Negative  Bias  (Volts) 

Figure  11:  Experimental  data  for  ground  and  flight  experiments 

for  each  correlated  area  independent  of  the  other  areas.  If  there  is  more  than  one  correlated 
area,  the  actual  arc  rate  for  the  array  is  the  sum  of  the  arc  rates  of  each  area. 

Cho  and  Hastings  use  this  procedure  in  Reference  3  to  calculate  the  arc  rate  numerically 
for  the  FIX  II  flight  and  ground  experiments.  As  can  be  seen  in  Figure  11,  the  results 
show  excellent  agreement  with  the  data  over  the  range  that  the  data  exists.  They  predict 
a  threshold  when  the  charging  process  is  exponentially  slow  and  also  predict  a  saturation 
for  high  voltages.  The  lower  parts  of  the  curves  cover  the  regime  where  the  enhanced 
field  electron  emission  charging  is  the  slowest  charging  process  in  the  system.  The  arc  rate 
dependence  on  voltage  here  is  exponential  and  enables  a  threshold  voltage  to  be  defined 
with  a  small  uncertainty.  This  threshold  voltage  can  be  defined  as  the  voltage  at  which  the 
arc  rate  is  decaying  very  rapidly.  The  upper  parts  of  the  arc  rate  curves  cover  the  regime 
dominated  by  the  ion  recharging  time.  This  leads  to  a  decrease  in  the  rate  of  change  of  arc- 
frequency  as  can  be  clearly  seen  in  the  data.  The  fact  that  the  arc  rate  scales  with  the  density 
for  the  higher  voltages  can  also  be  explained  from  the  dominance  of  the  ion  recharging  time 
since  this  scales  directly  with  density. 
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Figure  12:  Schematic  of  a  wrap-through-coniact  solar  array 


2.2  Wrap-Through-Contact  Cells 

A  schematic  of  a  wrap-through-contact  (WTC)  solar  cell  is  shown  in  Figure  12.  Kapton 
covers  the  metal  interconnector  so  it  is  not  exposed  to  the  ambient  plasma  environment 
like  the  interconnectors  of  conventional  cells.  One  of  the  reasons  for  this  design  was  to 
eliminate  arcing  at  the  interconnector-cell  interface.  On  the  edge  of  the  cell,  however,  the 
semiconductor  cell  itself  is  exposed.  Since  the  semiconductor  is  adjacent  to  both  a  dielectric 
coverglass  and  a  dielectric  substrate,  arcing  can  occur.  In  ground  tests,  arcing  occurred  on 
WTC  cells  at  bias  voltages  as  low  as  -400V.  Consequently,  a  model  is  needed  to  understand 
how  arcing  occurs  on  this  type  of  cell. 

As  with  the  conventional  cell,  the  area  of  interest  for  studying  electric  field  buildup 
can  be  simplified  to  two  dielectrics  and  a  conductor.  In  this  case,  the  conductor  is  situated 
between  the  two  dielectrics  as  shown  in  Figure  13.  The  numerical  model  for  the  conventional 
cells  could  be  modified  by  a  simple  change  of  boundary  conditions.  The  problem  is  more 
complex,  however,  as  the  conductor  is  now  in  the  computational  domain  instead  of  being 
merely  a  boundary  condition.  In  addition,  the  lower  dielectric  can  not  be  treated  as  a 
simple  boundary  as  the  conductor  was  in  the  conventional  cell  model.  To  properly  include 
the  dielectric  properties  of  the  substrate,  a  dielectric  of  two  grid  cells  thickness  is  added 
beneath  the  conductors  and  in  the  gap  between  the  cells.  The  new  geometry  also  has  two 
triple  Junctions  on  each  of  the  two  conductor  edges,  making  the  previous  grid  clustering 
inadequate.  The  grid  is,  therefore,  altered  to  again  cluster  near  the  triple  junctions  as  well 
as  along  the  side  surfaces,  as  shown  in  Figure  14. 

In  Section  2.1.1  the  results  of  the  ambient  ion  charging  calculated  by  the  PIC  scheme 


Coverglass 

(Dielectric) 


Si  Cell 
(Conductor) 


Triple  Junctions 


Kapton 

(Dielectric) 


Figure  13:  Wrap-through-contact  solar  array  model  used  for  numerical  simulations 


are  represented  by  the  electric  potential  plot  (Figure  6)  and  the  surface  charge  density  plot 
(Figure  7).  Corresponding  plots  for  the  WTC  geometry  are  shown  in  Figures  15  and  16. 
As  expected,  the  highest  potential  value  is  on  the  conductor  surface  with  gradients  falling 
off  quickly  around  this  voltage  source.  The  gap  is  sufficiently  large  that  the  gradients 
do  not  interfere  but  rather  connect  at  the  coverglass  surface.  The  potential  lines  fall  off 
uniformly  beyond  the  coverglass  surface.  Above  the  substrate  surface,  the  potential  gradients 
concentrate  near  the  conductor  with  potential  lines  peaking  sharply  in  the  center  due  to  the 
grid  configuration.  In  the  coverglass,  the  potential  gradients  are  curved  near  the  edge  of  the 
cell  but  straighten  away  from  the  cell  edge. 

To  simulate  the  electric  field  build  up,  the  modified  PIC  code  is  run  with  the  enhanced 
electron  field  emission  (EFEE)  charging  processes  included.  The  initial  conditions  are 
obtained  from  the  results  of  the  ambient  ion  charging  calculations  and  the  enhancement 
factor,  /?,  for  each  conductor  cell.  Due  to  the  high  electric  potential  normal  to  the  conductor 
surface,  much  lower  enhancement  factors  {(3  ^  30-60)  are  used  to  reduce  the  number  of 
electrons  emitted  from  the  conductor.  When  too  many  electrons  are  emitted,  the  EFEE 
charging  time  is  too  small,  making  it  less  than  or  on  the  order  of  a  capacitor  charging  time. 
The  timestep  also  affects  the  number  of  particles  injected  into  the  domain,  so  it  is  typically 
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Figure  15;  Typical  electric  potential  for  ambient  ion  charging  of  WTC  cells 


chosen  to  be  UpeAt  =  0.01.  Larger  timesteps  can  be  used  when  fewer  particles  are  in  the 
domain,  which  occurs  when  the  electric  field  at  the  triple  junction  has  not  increased  enough 
for  EFEE  charging  to  begin.  The  PIC  code  is  run  until  the  electric  field  runs  away  at  one 
or  both  of  the  triple  junctions.  Since  the  PIC  code  automatically  accounts  for  space  charge 
effects,  the  simulation  is  often  run  beyond  the  electric  field  runaway  into  the  space  charge 
current  limited  regime,  which  limits  the  electric  field  magnitude. 

The  cell  properties  used  for  the  WTC  simulations  are  based  on  the  Space  Station  Freedom 
WTC  cell.  The  coverglass  and  semiconductor  are  each  203^m  (8  mil)  thick,  and  the  cell 
gap  is  1mm.  The  semiconductor  is  silicon,  which  has  a  work  function  of  4.85eV.  The 
coverglass  is  assumed  to  be  ceria-doped  microsheet  (CMX)  with  a  dielectric  constant  of  4 
and  secondary  electron  properties  of  E,na.  -  400V  and  7r.«.x  =  4.  The  Kapton  substrate 
has  a  dielectric  constant  of  3.5  and  assumed  secondary  electron  properties  of  E„iax  -  300V 
3.nd  ’ymax  ~  3. 

The  EFEE  charging  of  the  WTC  cells  over  the  range  of  300-500V  is  distinguished  by  two 
classes  of  behavior,  the  first  occurring  with  lower  p  values  (~30)  and  the  second  occurring 
with  higher  /?  values  (~50-60).  As  seen  in  Figures  17  and  18,  the  electric  field  at  the  upper 
triple  Junction  increases  initially  for  the  first  class  and  decreases  initially  for  the  second 
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Figure  16:  Typical  surface  charge  density  along  the  side  dielecuic  surface  after  ambient  ion  charging  of  WTC 
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Figure  17:  Class  1  electric  field  at  upper  uiple  junction  versus  time 


class.  No  runaway  occurs  at  the  lower  triple  junction  within  this  time. 

In  the  first  class,  the  ambient  ion  charging  continues  to  build  up  the  electric  field  at  the 
upper  triple  junction  until  it  is  high  enough  to  initiate  EFEE  charging.  Once  initiated,  the  high 
flux  of  electrons  causes  the  field  to  decrease  for  a  short  time  before  the  runaway.  As  shown  in 
Figure  19,  the  surface  charge  along  the  side  of  the  coverglass  does  not  change  much  during 
the  ambient  ion  charging,  as  expected,  but  also  does  not  change  much  during  the  electric  field 
runaway.  The  surface  charge  along  the  front  surface,  however,  does  increase  substantially, 
indicating  that  the  electrons  from  the  conductor  are  striking  there  and  increasing  the  surface 
charge  through  secondary  electron  emission. 

In  the  second  class  of  behavior  EFEE  charging  begins  immediately,  emitting  many 
electrons  into  the  domain.  Although  the  surface  charge  density  does  increase,  as  shown  in 
Figure  20,  most  of  the  electrons  quickly  exit  the  domain  without  striking  any  of  the  cell 
surfaces.  Just  prior  to  runaway,  the  difference  between  the  number  of  electrons  emitted  and 
the  number  of  electrons  impacting  the  dielectric  increases  substantially  at  the  same  time 
that  the  total  number  of  electrons  in  the  domain  increases  substantially.  The  electric  field 
then  runs  away,  and  the  surface  charge  density  along  the  coverglass  side  and  front  surfaces 
increases  significantly. 
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Chapter  3 

Arc  Mitigation  Methods 


In  this  research  several  methods  of  reducing  the  arc  rate  were  studied.  As  expressed  in 
Eq.  (5),  the  arc  rate  is  the  inverse  of  the  sum  of  the  two  charging  times,  the  ambient  ion 
charging  time  and  the  enhanced  electron  field  emission  (EFEE)  charging  time  Tejee  ■ 
Since  non  is  dependent  mainly  on  mission  parameters  such  as  the  ambient  plasma  density,  arc 
rate  mitigation  can  only  be  achieved  by  increasing  Te/ec,  which  is  affected  by  cell  properties 
only.  From  Eq.  (7),  the  properties  which  affect  Te/ee  to  the  greatest  extent  are  those  in 
the  exponential  factor  and  the  secondary  electron  yield  which  causes  the  time  to  be 
non-existent  if  it  is  equal  to  or  less  than  unity.  In  addition,  lengthening  the  coverglass  over 
the  interconnector  to  obstruct  the  electron  trajectories  should  also  increase  Tefe.e,  if  only  by 
increasing  the  distance  over  which  the  surface  charge  must  build  up. 

For  all  numerical  simulations  described  in  this  chapter,  the  domain  size  used  is  3mm  in 
the  X  direction  and  2.5mm  in  the  y  direction,  which  includes  one-half  of  two  2mm  wide 
cells  and  a  0.5mm  gap  in  between  them,  as  shown  in  Figure  5.  The  simulations  assumed  the 
same  environment  of  n,,  =5xl0''m-3,  T,  =  T,.  =0.1eV,  and  a  kinetic  energy  of  incoming 
ions  of  5eV,  all  typical  of  low  earth  orbit.  A  90‘’  orientation  to  the  ram  velocity  was  used 
for  simphcity.  In  addition,  the  emission  site  on  the  conductor  was  set  adjacent  to  the  triple 
junction  with  an  area  of  Sreai  =E2x  10" as  determined  by  the  grid  cell  length  at  that 
location.  Since  EFEE  charging  is  dependent  on  an  emission  site  near  the  triple  junction, 
this  condition  should  define  the  upper  bound  for 

3.1  Control  Case 

The  cell  used  as  the  control  case  for  these  arc  rate  simulations  is  the  silicon  conventional 
cell  without  a  coverglass  overhang.  The  input  parameters  chosen  to  simulate  this  cell  are 
shown  in  Table  1.  The  dielectrics  are  a  fused  silica  coverglass  and  DC  93500  adhesive,  and 
the  interconnector  material  is  Kovar.  The  interconnector  work  function  is  taken  to  be  the 
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Table  1:  Conventional  silicon  cell  data  used  in  numerical  simulations 


d, 

0.153mm 

0.037mm 

3.5 

2.7 

')max\ 

3.46 

'^inax2 

3.0 

EftiaxX 

330V 

Efjiax2 

300V 

dgo/p 

0.5mm 

4>w 

4.76eV 

weighted  average  of  the  work  functions  of  the  elements  which  compose  it.  The  dielectric 
thicknesses  are  typical  values. 

To  determine  the  effect  of  varying  each  parameter,  two  calculations  are  made.  First,  the 
numerical  code  discussed  in  Section  2. 1 . 1  is  used  to  determine  Tc  fee  over  a  range  of  typically 
5  values  of  dV.  The  charging  time  Te/ce  is  then  plotted  against  i3V  since  the  analytic  theory^ 
indicates  that  Tefee  is  a  strong  function  of  PV.  This  plot  is  shown  in  Figure  21  for  the  control 
case  over  a  relevant  range.  Times  less  than  1  x  10"®sec  are  not  useful  as  that  interval  is  on 
the  order  of  a  capacitor  build-up  time.  Times  greater  than  1  x  lO'^sec  are  also  not  useful  as 
either  the  ion  charging  time  would  dominate  or  the  orbit  would  be  completed,  causing  the 
power  system  to  regenerate.  In  the  following  property  variations,  however,  the  same  range 
of  PV  is  maintained  where  possible  to  simplify  comparisons. 

Second,  the  analytical  model  discussed  in  Section  2.1.2  is  used  to  determine  the  arc 
rates.  In  this  model,  the  arcing  time  is  determined  by  the  one  emission  site  on  a  conductor 
of  typically  1000  sites  that  has  the  shortest  charging  time  of  all  the  cells.  This  is  a  more 
realistic  simulation  than  the  numerical  simulations  in  which  only  the  effect  of  one  site  on  one 
conductor  was  studied.  Hence,  the  effects  determined  by  the  numerical  simulations  often 
do  not  have  as  much  impact  on  the  arc  rate  results  as  one  might  expect.  The  experiment  time 
is  also  a  consideration,  particularly  for  lower  voltages  where  arcing  takes  a  longer  time  to 
occur.  For  the  control  case  and  the  mitigation  cases  the  experiment  time  is  arbitrarily  chosen 
to  be  one  second.  Also,  arc  rates  are  calculated  at  intervals  of  -  lOOV.  The  arcing  rates  for  the 
control  case  are  shown  in  Figure  22.  The  curve  clearly  shows  the  two  dominating  regions  of 
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Tejee  and  At  lower  voltages,  T^fcc  dominates  and  the  arc  rate  is  a  strong  function  of  the 
voltage.  At  higher  voltages,  Tian  limits  the  arcing  rate  to  the  value  determined  by  mission 
parameters,  which  is  a  weak  function  of  the  voltage.  Consequently,  the  arc  rate  mitigation 
methods  studied  are  intended  to  shift  and  alter  the  slope  of  the  Te/ee  dominated  region.  The 
Tton  dominated  region  can  also  be  shifted  though  it  will  remain  at  the  same  arc  rate.  The  arc 
rate  results  are  presented  at  the  end  of  the  chapter  in  a  separate  section  so  that  comparisons 
may  be  made  among  all  methods  studied. 


3.2  Interconnector  Material 


The  work  function  of  the  electron  emitting  surface  determines  the  ease  with  which  electrons 
are  released.  If  the  number  of  electrons  emitted  from  the  interconnector  is  reduced,  the  time 
for  the  electric  field  at  the  triple  junction  to  build  up  will  be  increased.  In  the  analytical  for¬ 
mula  for  Te/ee,  the  work  function  determines  the  value  of  the  Fowler-Nordheim  coefficients 
A  and  B  given  by  Eq.  (8)  and  (9).  Figure  23  shows  Te/ee  plotted  against  dV  over  a  range 
of  significant  values  for  varying  work  functions.  As  expected  metals  with  work  functions 
higher  than  the  control  case  of  4.76eV  have  longer  times  for  EFEE  charging,  and  metals 
with  lower  work  functions  have  shorter  times. 

These  numerical  simulation  results  can  easily  be  predicted  from  the  theory.  To  determine 
the  effect  of  a  different  work  function  we  can  solve  the  ratio  t^, fee{<Pw2  V  TejetiPw )  using 

Eq.  (7); 

’^e/ce(0tti2 )  ( X  10  ci 

=  exp  (0^2  -  (Pw  ) 


Tefeci^Pu,)  ^  J 

The  only  unknown  variable  is  r/,  which  is  within  the  range  1.001-1.005  for  emission  sites 
adjacent  to  the  triple  junction.  For  the  cases  shown  in  Figure  23,  this  expression  is  evaluated 
and  plotted  with  the  equivalent  numerical  values  in  Figure  24.  The  predictions  are  all  within 
the  margin  of  error  for  the  analytical  and  numerical  values,  except  for  the  extreme  case  of 
=5.9eV  which  is  predicted  about  an  order  of  magnitude  too  high.  This  is  therefore  a 
useful  tool  for  predicting  the  effect  of  changing  the  exposed  interconnector  metal. 


3.3  Dielectric  Thickness 

Another  method  of  reducing  the  arc  rate  is  to  increase  the  dielectric  thickness,  d  =  +  d2. 

This  increases  the  surface  on  which  charge  can  accumulate  and  reduces  the  average  charge 
density.  The  analytic  theory®  accounts  for  changes  in  thickness  in  the  exponential  factor  of 
Tefce-  The  exponential  factor  also  includes  which  is  dependent  on  the  ratio  of  dielectric 


33 


pv  X  10®  (volts) 

Figure  23:  Enhanced  field  electron  emission  charging  time,  Tg/ee-  versus  0V  for  different  work  functions, 
(eV) 


pv  X  10®  (volts) 


Figure  24:  Analytic  predictions  and  numerical  results  for  Tf~ffc/T(-fee{(l>w  =  4.76eK )  versus  0V 


Figure  25:  Enhanced  field  electron  emission  charging  lime,  Tf  /er,  versus  0V  for  different  dielectric  thick¬ 
nesses,  diixm) 

thicknesses,  di/d2.  In  order  to  exclude  this  factor,  the  ratio  is  the  same  for  each  case 
studied.  Further  simulations  which  studied  the  effect  of  this  ratio  showed  that  it  did  not  alter 
the  results.  Figure  25  shows  Tefee  plotted  against  0V  for  the  control  case  and  for  four  cases 
with  differing  thicknesses.  As  expected,  increases  in  thickness  increase  Tcfee  significantly. 

By  calculating  the  ratio  of  r^ff  Ad  -t-  Ad) /Tcfc^id)  using  Eq.  (7),  the  change  in  T(,f tf  can 
again  be  easily  predicted: 


Tcfec{d+Ad) 

'  f  d  +  Ad 

j  —  exp 

[[  d  ) 

\0V7]ioJ\ 

The  analytic  results  are  plotted  with  the  numerical  results  in  Figure  26.  The  predictions  are 
all  within  the  accuracy  of  the  numerical  simulations.  The  prediction  for  the  largest  variation 
of  ci  =  118//m  is  only  off  by  a  factor  of  3.  Using  the  same  control  case,  an  increase  in 
thickness  to  250//m  would  result  in  an  increase  in  Tf/cf  of  nearly  6  orders  of  magnitude. 
This  will  effectively  eliminate  electric  field  run  away.  Meanwhile,  a  decrease  in  thickness 
to  50/.im  would  result  in  a  decrease  in  Tcfee  of  nearly  14  orders  of  magnitude.  In  this  case 
the  arcing  rate  will  be  dominated  by  the  ion  recharging  time,  ^0,*.  These  results  indicate 
that  the  new  thinner  solar  cells  such  as  Advanced  Photovoltaic  Solar  Array  (APSA)  will 
experience  higher  arcing  rates  at  relatively  low  voltages. 
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Figure  26:  Analytic  predictions  and  numerical  results  for  Tcff.el'Tcfcdd  =  190/U7n)  versus  0V 

3.4  Secondary  Electron  Yield 


The  secondary  electron  yield,  from  the  dielectrics  must  be  greater  than  unity  for  the 
electric  field  to  build  up  and  run  away.  If  is  equal  to  unity,  the  electric  field  will  assume 
a  steady  state  and  no  charge  accumulation  will  occur.  If  7^^  is  less  than  unity,  the  electric 
field  will  decrease  as  negative  charge  accumulates  on  the  dielectrics.  This  is  modeled  in 
Eq.  (7)  by  the  factor  (7sf  -  l)~'. 

A  sensitivity  scan  of  over  relevant  values  of  dV  for  7«  values  of  actual  materials  is 
shown  in  Figure  27.  As  in  the  previous  sensitivity  scans,  only  the  property  of  interest,  7„  , 
is  varied  from  the  control  case.  In  these  cases,  7,,.  is  the  same  for  both  dielectrics,  although 
the  control  case  has  slightly  different  values  of  7ef  for  the  two  dielectric  materials,  fused 
silica  and  DC  93500.  In  order  to  choose  relevant  values,  the  Eq.  (34)  for  7^  is  evaluated 
for  different  7,^^  and  values  determined  for  actual  materials.  Since  the  ambient  ion 
charging  is  not  calculated  with  secondary  electron  effects,  the  initial  conditions  for  EFEE 
charging  are  the  same  as  for  the  control  case  for  every  7^^  variation.  Consequently,  d, 
must  only  be  known  for  the  control  case  to  determine  7„  for  any  other  set  of  secondary 
electron  parameters.  From  the  numerical  simulation  of  the  control  case,  rf,  is  35.6/im,  so 
9^  is  estimated  to  be  89'’,  and  E^  is  estimated  to  be  75V.  For  the  control  case  the  estimated 
7ef  is  15,  but  the  numerical  results  showed  7ef  to  be  7.  In  the  other  cases  studied  the 
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Figure  27:  Enhanced  field  electron  emission  charging  time,  Tf  /ee,  versus  0V  for  different  secondary  electron 
yields,  ')ee. 

estimated  values  of  7e£  are  also  about  a  factor  of  2  higher  than  the  calculated  values.  The 
lowest  estimated  value  studied  is  =  2.5  for  cesium’s  properties,  but  the  calculated  value 
is  7ee  =  0.7- 1.2.  The  latter  corresponds  to  the  result  of  no  electric  field  run  away  by  the 
numerical  simulation.  The  results  are  shown  in  Figure  27  with  the  numerical  values  of 
7ee.  As  expected,  cases  with  secondary  electron  yields  near  unity  have  much  longer  EFEE 
charging  times,  and  the  case  (Cs)  with  7^0  values  near  and  less  than  unity  has  an  infinite 
EFEE  charging  time  since  the  electric  field  at  the  triple  junction  did  not  run  away. 

One  difference  between  the  analytical  and  numerical  results  is  the  effect  of  the  dielectric 
surface  charge,  which  is  not  taken  into  account  in  the  analytical  model.  As  shown  in 
Figure  28,  the  surface  charge  next  to  the  triple  junction  is  negative  during  EFEE  charging 
for  the  Cs  case  of  7ef  =  0.7- 1.2,  causing  the  electric  field  to  be  unable  to  run  away.  In  the 
Li  case  of  7ef.  =  0.4-2,  the  surface  charge  builds  up  over  time  as  it  must  for  the  electric  field 
to  build  up.  This  surface  charge  density  effect  is  discussed  further  in  the  next  section. 

3.5  Dielectric  Constants 

The  effect  of  different  dielectric  constants  is  not  as  apparent  in  the  analytical  formula  for 
Tefee  (Eq.  (7))  as  for  the  other  cell  properties.  In  the  analytical  formula,  the  dielectric 
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Figure  29:  Enhanced  field  electron  emission  charging  time,  t^, fgf  ,  versus  0V  for  different  dielectric  constant 
ratios 

constants  only  directly  affect  ^  and  ^o-  The  ratio  in  Eq.  (10)  determines  whether 

^0  is  less  than,  equal  to,  or  greater  than  1 .  Although  this  does  not  significantly  affect  the 
exponential  factor  of  Tg/ee,  it  does  affect  7fc  through  E^  as  shown  in  Eq.  (34). 

The  dielectric  constants,  however,  significantly  affect  the  surface  charge  density  and 
consequently  the  secondary  electron  yield.  Since  td,Ei  =  td^E2,  the  ratio  of  dielectric 
constants,  tdjtd^,  must  equal  the  ratio  of  electric  fields,  E-z/E-^.  Therefore,  if  tdjtd^  is 
greater  than  one,  the  electric  field  on  the  lower  side  surface  will  be  higher.  This  higher 
electric  field  results  in  a  higher  surface  charge  density  causing  more  electrons  to  be  attracted 
and  many  more  secondary  electrons  to  be  emitted.  The  electric  field  at  the  triple  junction 
should  then  build  up  quite  rapidly  causing  r^fcc  to  be  relatively  short.  If  td,  /td^  is  less  than 
one,  the  electric  field  on  the  lower  surface  will  be  lower  than  the  electric  field  on  the  upper 
surface,  causing  it  to  be  very  difficult  for  the  electric  field  to  build  up  at  the  triple  Junction. 
The  EFEE  charging  time  should  then  be  very  long. 

The  effect  of  the  surface  charge  density  is  not  accounted  for,  however,  as  the  analytic 
expressions  do  not  predict  the  large  variation  of  TcJcc  calculated  by  the  numerical  simulations. 
The  results  are  shown  in  Figure  29.  By  increasing  the  ratio  td, /td^  from  the  control  case  of 
1.3  to  2.7,  Te/ff  decreases  by  about  7  orders  of  magnitude.  By  decreasing  the  ratio  from  1.3 
to  0.74,  Te/ee  increases  by  about  12  orders  of  magnitude.  As  expected,  these  results  can  be 
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distance  from  triple  junction  (mm) 

I  Figure  30:  Dielectric  side  surface  charge  density  before  EFEE  charging  for  different  dielecu-ic  constants 

i 

j 

explained  by  the  surface  charge  density.  The  numerical  results  of  the  surface  charge  density 
calculations  are  shown  in  Figure  30.  For  higher  values  of  e^,  /e^j ,  the  surface  charge  density 
is  correspondingly  high  over  the  lower  dielectric  side  surface.  This  attracts  electrons,  as 
shown  in  Figure  31,  causing  the  electric  field  to  run  away  in  a  short  time.  For  the  control 
case,  the  surface  charge  density  does  not  build  up  as  high  before  EFEE  charging  so  T(Jc( 
is  longer.  Correspondingly,  the  electron  trajectories  are  not  as  concentrated  near  the  triple 
junction,  as  shown  in  Figure  32.  For  cases  with  tdjtd^  <  1,  the  surface  charge  density 
over  the  lower  dielectric  side  surface  is  negative  before  beginning  EFEE  charging  so  that 
electrons  are  repelled  from  the  adhesive,  as  shown  in  Figure  33,  causing  secondary  electron 
emission  to  be  substantially  reduced.  For  these  latter  cases  the  useful  secondary  electron 
yield  is  near  or  below  1 ,  so  that  the  electric  field  assumes  a  steady  state  which  does  not  run 
away  in  any  reasonable  time. 
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Figure  32:  Electron  u-ajectories  for  the  control  case  of  Cd,  j td^  —  \  .3 
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Figure  33:  Electron  trajectories  for  tdjtd^  —  0.74 

The  effects  on  the  secondary  electron  yield  as  calculated  by  the  numerical  simulations 
can  be  seen  in  Figure  34.  As  expected  from  the  surface  charge  density  results,  7ee  for  cases 
with  td,  <  1  is  significantly  lower  than  for  the  other  cases. 

3.6  Overhanging  the  Coverglass 

The  final  mitigation  strategy  studied  involves  lengthening  the  coverglass  into  an  overhang 
over  the  interconnector  to  create  a  back  surface,  as  shown  in  Figure  35,  on  which  charge 
can  accumulate.  This  method  was  expected  to  increase  since  the  overhang  creates  an 
additional  surface  over  which  the  charge  must  build  up  for  the  electric  field  to  run  away. 

3.6.1  Numerical  Results 

As  shown  in  Figure  36,  the  results  confirm  this  expectation.  As  the  overhang  is  increased  in 
length  up  to  30nm,  T^fee  is  increased  uniformly.  In  the  range  studied,  this  pattern  changes 
for  overhangs  longer  than  30//m,  indicating  a  different  physical  effect  of  the  overhang  on 
the  EFEE  charging  time. 

The  difference  between  these  two  classes  is  clearly  shown  in  Figure  37.  The  electric 
field  at  the  triple  junction  Etj  increases  with  time  for  the  the  first  class,  indicating  a  build 
up  of  charge  due  to  EFEE  and  secondary  electron  emission.  However,  it  decreases  initially 
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Figure  34:  Average  secondary  electron  yield  over  the  adhesive  versus  time  for  different  dielectric  constants 
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Figure  35:  Model  of  coverglass  overhang 
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Figure  36:  Enhanced  field  electron  emission  charging  time,  Tg/ee,  versus  (5V  for  different  overhang  lengths 

before  increasing,  often  drastically,  for  the  second  class,  indicating  an  inital  build  up  of 
negative  charge  from  electron  accumulation  before  electron  emission  from  the  side  surface 
becomes  high  enough  to  create  the  positive  charge  build  up  necessary  for  runaway. 

The  corresponding  effects  of  these  differing  electric  fields  on  the  electron  trajectories 
can  be  seen  in  Figures  38  and  39.  The  first  shows  that  the  electrons  in  the  first  class  of 
simulations  are  merely  diverted  by  the  overhang.  The  second  shows  that  the  electrons  in  the 
second  class  are  mainly  confined  to  the  back  surface  with  few  arriving  at  the  side  surface  of 
the  coverglass.  This  shows  that  the  first  class  of  overhangs  are  acting  as  additional  surface 
area  while  the  second  class  are  obstructing  the  electron  trajectories  to  the  side  surface. 

The  potential  along  the  surface  reflects  these  electron  trajectories.  Figure  40  shows  that 
for  the  first  class  the  potential  along  the  back  surface  increases  with  time,  which  causes 
the  electric  field  at  the  triple  junction  to  increase  with  time  as  well.  For  the  second  class, 
shown  in  Figure  4 1 ,  the  back  surface  potential  initially  decreases  with  time  until  it  eventually 
increases  at  very  large  times.  This  initial  potential  decrease  corresponds  to  the  initial  elecu-ic 
field  decrease,  indicating  that  the  surface  is  charging  negatively  due  to  electron  build-up 
rather  than  positively  due  to  high  secondary  electron  emission.  This  is  clearly  the  case  as 
shown  in  Figure  39,  in  which  almost  all  of  the  electrons  strike  the  back  surface  while  very 
few  reach  the  side  surface.  Hence,  the  large  negative  charge  on  and  near  the  back  surface 


44 


Electric  Field  al  Triple  Junction  (V/m) 


29.704  29.836  29.692  29.946  30.000 


X  (mm) 

Figure  39:  Electron  trajectories  for  50//in  overhang;  0V  =  3x  lO^V 
suppresses  the  electric  field,  causing  it  to  decrease  rather  than  increase. 


3.6.2  Analysis 


The  results  above  suggest  that  there  is  a  critical  overhang  which  determines  the  behavior  of 
the  solution.  Below  this  critical  overhang,  the  electric  field  runaway  is  similar  to  the  case 
of  no  overhang  but  at  a  lengthened  time.  Above  it,  the  field  initially  decreases  for  a  long 
time  before  building  up  to  the  runaway.  Hence,  coverglasses  made  longer  than  this  critical 
overhang  will  have  substantially  reduced  arcing. 

Using  the  geometry  shown  in  Figure  2,  a  simple  calculation  can  determine  the  approx¬ 
imate  critical  overhang.  Since  the  electric  field  is  only  orthogonal  to  the  conductor  and 
approximately  orthogonal  to  the  back  surface,  the  free  electron  force  balance  is 

dvx  _ 

me—j^=  eEr  (44) 

and 

dvy 

where  me  is  the  electron  mass,  Vx  and  Vy  are  the  electron  velocities  in  the  x  and  y  directions 
respectively,  and  Ex  is  the  electric  field  in  the  x  direction. 

Electrons  emitted  from  the  back  surface  follow  ballistic  trajectories  according  to  this 
electric  field.  The  time  to  the  peak  of  each  ballistic  trajectory  is 

Vn 


t  = 


■'Xq 


eExInie ' 


(4(i) 
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where  is  the  initial  electron  velocity  in  the  x  direction.  The  distance  the  electron  travels 
along  the  back  face  before  re-impacting  the  surface  is 


yt 


(47) 


eEj^/nie  ’ 

where  Vy^  is  the  initial  electron  velocity  in  the  y  direction. 

The  electric  field  Ex  can  be  found  from  the  potential  difference  between  the  conductor 
and  the  back  surface.  Assuming  the  potential  of  the  front  surface  is  zero,  the  electric  fields 
of  the  coverglass  and  adhesive  can  be  expressed  as 


(48) 


and 


(49) 


where  (f)i  is  the  potential  along  the  interface  of  the  coverglass  and  adhesive  and  is  the  po¬ 
tential  of  the  conductor,  or  interconnector.  At  the  interface,  td,  Ex  ^  td^E^.  By  substitution 
to  eliminate  0,  and  E^ ,  both  of  which  are  unknown,  the  following  expression  is  derived  for 
E2,  or  Ex'. 


Ex  =  ^ 


(■d,d2l  td^dx 


(50) 


1  +  e<ijd2/ j 

Since  the  critical  overhang  determines  the  behavior  of  the  system,  the  differences  between 
the  two  classes  should  be  reviewed.  For  the  first  class,  the  electric  field  at  the  triple  junction 


increases  with  time.  This  means  that  positive  charge  is  accumulating  on  the  back  surface  of 
the  coverglass,  as  shown  in  Figure  40.  Therefore  the  secondary  electrons  are  escaping  from 
the  back  surface,  as  shown  in  Figure  38.  In  contrast,  in  the  second  class  of  solutions  the 
electric  field  is  initially  decreasing  with  time;  therefore,  negative  charge  is  accumulating  on 
the  back  surface  of  the  coverglass,  as  shown  in  Figure  41,  so  secondary  electrons  cannot 


immediately  escape  the  back  surface,  as  shown  in  Figure  39. 

Hence,  the  critical  overhang  that  separates  the  two  solutions  can  be  bounded  by  requiring 
that  the  distance  yt  <  do  so  that  electrons  cannot  immediately  escape  and  by  requiring  that 
the  energy  ruel^Vy^Vx^  be  bounded  by  the  secondary  electron  emission  energy  for  unity 
yield.  Ex-  This  latter  requirement  means  that  an  electron  striking  the  back  surface  never  has 
enough  energy  to  release  more  than  one  electron.  Therefore,  the  back  surface  must  charge 
negatively.  These  two  bounds  give 


dp  ^  4^1  1  4-  tdtd2ltd^dx 

di  ~  e0c  td,d2ltd^di  ■ 

For  the  conditions  in  Table  \,Ex^  40eV  and  with  0^  =  400V,  Eq.  (51)  gives  d^  <  61  ym. 
This  corresponds  well  with  what  is  observed  in  Figure  36  where  the  50/im  overhang  shows 
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the  characteristics  of  both  classes  while  the  75/rin  overhang  shows  only  the  characteristic 
of  the  second  class. 

This  approximate  expression  can  be  further  simplified  in  the  limit  ^2  <  d\  (which  is 
usually  the  case).  Eq.  (51)  then  simplies  to 


jc  ^  I 

do  <  —7 - fll- 

e<f>c  Crfi 


(52) 


This  shows  that  one  can  obtain  the  smallest  critical  overhang  by  reducing  d^  as  much 
as  possible,  consistent  of  course  with  the  need  to  maintain  radiation  protection,  and  by 
modifying  the  secondary  electron  emission  properties  of  the  coverglass  so  as  to  reduce  the 
energy  at  which  the  yield  is  unity. 


3.7  Arc  Rate  Results 


The  arc  rates  are  determined  by  using  the  analytical  model  discussed  in  Section  2.1.2,  as 
explained  for  the  control  case.  This  model  is  limited  by  the  analytic  expressions,  which  are 
only  approximations  of  real  processes  and  effects.  For  the  results  in  the  previous  section 
which  are  consistent  with  the  analytic  predictions,  the  model  is  used  with  the  varying 
parameters.  For  the  secondary  electron  yield,  dielectric  constant,  and  overhang  effects  on 
Tf/ee,  which  are  not  fully  accounted  for  in  the  analytic  expression,  the  model  is  used  with 
Te/ee  modified  according  to  the  results  from  the  numerical  simulations. 

The  results  are  shown  in  Figures  42-46.  In  all  cases  the  curves  saturate  around  350sec“ ' , 
where  the  arc  rate  is  hmited  by  Tion  which  is  not  as  dependent  on  voltage  as  is  Te/ee-  The 
goal  of  the  mitigation  strategies  can  clearly  be  seen,  then,  as  increasing  the  voltage  range  of 
Te fee  dominance. 

The  effect  of  different  work  functions  on  the  arc  rate  is  shown  in  Figure  42.  None  of  the 
work  functions  significantly  affect  the  arc  rate,  although  the  higher  work  functions  do  not 
arc  at  -2()0V.  Also,  little  difference  is  seen  between  the  arc  rates  of  commonly  used  Kovar 
and  silver  interconnectors. 

The  thickness  variations  show  a  greater  effect  on  the  arc  rate,  especially  for  the  smaller 
thicknesses  between  -200V  to  -400V.  For  a  1 18^m  thickness  of  the  coverglass  and  adhesive 
dielectrics,  the  arc  rate  is  almost  saturated  at  -300V,  while  the  control  case  thickness  of 
190/im  does  not  saturate  until  -600V.  Increasing  the  thickness  to  200//m  or  even  226/xm 
eliminates  arcing  at  -200V  but  does  not  reduce  arcing  significantly  at  higher  voltages. 

The  effect  of  the  secondary  electron  yield  on  the  arc  rate  is  shown  to  be  significant  only 
when  the  electric  field  does  not  run  away,  as  in  the  jee  =  0.7- 1.2  case.  Even  for  the  two 
cases  (BeO  and  Li)  which  have  secondary  electron  yields  near  and  below  unity,  the  arc  rate 
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Figure  42:  Analytic  arc  rates  for  varying  interconnector  work  functions,  (pw 


Vott^e 


Figure  43:  Analytic  arc  rates  for  varying  dielectric  thicknesses,  d 
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is  just  barely  reduced  from  the  control  case.  Any  secondary  electron  yield  not  near  unity, 
such  as  that  of  carbon  (%e.  =  4-4.5),  did  not  change  the  arcing  rates  at  all  from  the  control 
case  of  jec  =  7. 

The  dielectric  constant  variations  and  coverglass  overhang  cases  show  the  largest  effect 
on  the  arc  rate.  For  the  =  2.7  case,  the  arc  rate  is  almost  saturated  at  -300V,  For  the 
case  of  equal  dielectric  constants,  however,  the  arc  rate  does  not  saturate  until  -900V.  For 
variations  of  td^/td^  <  1  the  arc  rate  is  reduced  significantly.  No  arcing  is  predicted  until 
-600V  for  the  td,/td2  =  0.74  case. 

For  the  coverglass  overhang  arc  rates,  the  two  different  classes  can  be  seen  clearly  in 
Figure  46.  The  arc  rates  of  the  first  class  do  not  decrease  much  even  with  an  overhang  of 
30/j,m,  but  the  arc  rates  of  the  second  class  are  significantly  reduced.  The  75/im  case,  which 
is  beyond  the  critical  overhang  for  the  cell  studied,  shows  no  arcing  until  -600V  and  does 
not  reach  the  saturation  level  in  the  range  of  voltages  examined. 
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Figure  46:  Predicted  arc  rates  for  varying  coverglass  overhang  lengths,  dp  (/rm) 
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Chapter  4 

PASP  Plus  Predictions 


In  March,  1994,  the  Solar  Array  Module  Plasma  Interactions  Experiment  (SAMPIE)  was 
flown  aboard  the  Space  Shuttle.  In  the  next  few  years,  three  more  space  missions  are 
planned  to  study  high  voltage  solar  arrays,  including  the  Photovoltaic  Array  Space  Power 
Plus  Diagnostics  (PASP  Plus)  experiment.  By  using  the  numerical  and  analytical  models 
for  the  conventional  and  wrap-through-contact  cells,  arcing  predictions  were  made  for  the 
PASP  Plus  mission  to  aid  in  it’s  preparation  and  to  further  check  the  vahdity  of  these  models. 


4.1  PASP  Plus 

PASP  Plus  is  designed  to  provide  the  Air  Force  with  important  information  regarding  higher 
power  space  systems.  Since  greater  spacecraft  power  will  be  required,  the  Air  Force  an¬ 
ticipates  using  advanced  high  voltage  solar  arrays  and  power  distribution  systems.  The 
environmental  interaction  hazards  of  these  systems  therefore  need  to  be  determined  to  im¬ 
prove  the  reliability  of  future  missions.  To  accomplish  its  task,  the  PASP  Plus  experiment 
includes  both  advanced  solar  arrays  and  a  number  of  diagnostic  instruments.  PASP  Plus 
has  four  main  objectives^®: 

(1)  characterize  the  electrical  performance  and  the  environmental  interaction  of  advanced 
solar  array  designs  operating  at  high  voltages  in  the  ionosphere  and  magnetosphere; 

(2)  determine  the  long-term  radiation  degradation  effects  of  several  advanced  solar  arrays; 

(3)  determine  the  impact  of  the  space  environment  on  solar  array  operation  for  various 
solar  cell  technologies;  and 

(4)  obtain  flight  performance  data  for  advanced  array  designs  never  before  flown. 
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4.2  Experiment  Description 


The  PASP  Plus  experiment  will  be  flown  on  the  Advanced  Photovoltaic  and  Electronics 
Experiments  (APEX)  satellite,  scheduled  to  be  launched  in  the  summer  of  1994  aboard 
the  Pegasus  launch  vehicle.  The  satellite  should  reach  an  orbit  of  1950km  apogee,  360km 
perigee,  and  70"  inclination.  Throughout  the  orbit,  the  spacecraft  will  to  be  continuously 
sun  pointing.  The  experiment  is  expected  to  last  one  year  minimum  and  three  years  nominal. 

The  APEX  spacecraft  in  its  deployed  configuration  is  shown  in  Figure  47.  The  arrays 
are  located  on  the  top  shelf  and  on  one  of  the  deployed  panels.  The  other  three  panels  are 
the  solar  arrays  needed  for  the  spacecraft  power  system.  Beneath  the  payload  shelf  is  the 
avionics  shelf,  where  most  of  the  diagnostic  instruments  are  located.  There  are  three  types 
of  sensors  to  obtain  array  performance  data  and  five  types  of  sensors  to  characterize  the 
space  environment  around  the  spacecraft.  To  calibrate  the  performance  of  the  solar  arrays, 
the  following  instruments  will  be  used:  a  sun  sensor  to  measure  the  incidence  angle,  meters 
to  measure  the  current  and  voltage,  and  sensors  near  each  array  to  measure  the  temperature. 
Instruments  which  characterize  the  space  environment  include  a  transient  pulse  monitor  for 
measuring  electromagnetic  interference,  a  langmuir  probe  with  sense  potential  capability  for 
measuring  the  thermal  plasma  temperature  and  density,  a  dosimeter  for  measuring  ion  and 
electron  radiation,  an  electrostatic  analyzer  to  determine  electron  and  ion  spectra  of  moderate 
energy  (lOeV  -  30keV),  and  five  monitors  for  determining  the  amount  of  degradation  due 
to  radiation  and  the  amount  due  to  contamination.  The  spacecraft  also  contains  a  voltage 
generator  capable  of  biasing  the  arrays  at  multiple  high  voltage  levels,  an  electron  emitter 
which  will  enable  the  satellite’s  potential  to  be  altered,  and  an  electronics  controller  for 
satellite  operation. 

Seventeen  arrays  of  twelve  different  types  will  be  mounted  on  the  payload  shelf  and  the 
deployed  panel.  Since  there  are  only  16  electrical  channels,  one  of  the  arrays  will  not  be 
used.  Ten  of  these  arrays  will  be  subject  to  multiple  voltage  biasing.  The  arrays,  as  well  as 
the  instruments  described  above,  are  all  numbered  for  identification.  They  will  be  therefore 
introduced  in  order  of  their  instrument  numbers,  and  the  biased  arrays  will  then  be  described 
in  further  detail,  with  the  exception  of  the  concentrator  arrays  which  are  not  included  in 
this  research.  Arrays  #0,  #1,  and  #2  contain  silicon  conventional  cells.  The  silicon  WTC 
cells  designed  for  Space  Station  Freedom  make  up  array  #3.  Arrays  #4  -  #6  are  Applied 
Solar  Energy  Corporation’s  (ASEC)  gallium  arsenide/germanium  (GaAs/Ge)  conventional 
cells,  but  array  #5  will  not  be  used  so  that  its  channel  can  be  used  by  the  APSA  ceU  array. 
The  advanced  solar  cells  of  AlGaAs/GaAs  Monolithic  MBG  are  on  array  #7.  Spectrolab’s 
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Figure  47:  Deployed  APEX  spacecraft  with  PASP  Plus  experiment  payload 


WTC  GaAs/Ge  cells  are  located  on  array  #8.  Array  #9  contains  the  amorphous  silicon  cells 
developed  by  TRW  and  Solarex.  Another  advanced  solar  cell,  indium  phosphate  (InP), 
is  on  array  #10.  The  conventional  GaAs/Ge  cells  made  by  Spectrolab  are  on  array  #1 1. 
Another  MBG  cell,  Boeing’s  GaAs/CuInSe  solar  cell,  make  up  arrays  #12  and  #13.  The 
Mini-Cassegrainian  GaAs  cell  created  at  TRW  will  be  on  array  #14  to  accompany  the  newer 
Mini-Dome  Fresnel  Concentrator  created  at  Boeing  which  is  on  array  #15.  Finally,  the 
APSA  cell  was  added  late  to  the  list  so  its  array  is  #36. 

The  silicon  conventional  arrays  #1  and  #2  are  the  first  biased  arrays,  containing  20 
and  60  cells  respectively.  Originally  built  by  RCA  in  1984,  these  cells  are  included  to  be 
representative  of  present  operational  space  flight  solar  cells.  Also,  they  were  not  designed 
to  withstand  high  voltage  operation.  The  cells  are  2cmx4cm  and  203/.im  (8mil)  thick.  The 
coverglass  is  fused  silica  and  the  cover  adhesive  is  DC  93500,  of  estimated  thickness  1 53/im 
(3mil)  and  37^m  (1.5mil)  respectively.  The  interconnector  is  assumed  to  be  Kovar,  which 
has  an  estimated  work  function  4.76eV. 

The  Space  Station  Freedom  array  of  4  cells  (#3)  will  also  be  biased.  These  WTC  cells 
were  designed  to  operate  in  low  earth  orbit  at  a  160V  nominal  operating  voltage.  The 
8cm  X  8cm  cells  are  203//m  (8mil)  thick.  As  discussed  in  Section  2.2,  the  cells  are  covered 


with  a  ceria-doped  borosilicate  coverglass  and  a  Kapton  substrate. 

Arrays  #4  and  #6  of  ASEC’s  GaAs/Ge  conventional  cell  design  are  comprised  of  20 
and  12  cells  respectively,  and  both  will  be  biased.  These  cells  were  designed  to  have  high 
radiation  resistance  and  high  cell  efficiency,  making  them  representative  of  the  standard  for 
the  next  generation  of  satellite  photovoltaic  power.  The  4cmx4cm  cells  have  a  thickness  of 
76/.im  (3.5mil)  and  CMX  coverglasses  of  102/im  (4mil).  The  interconnector  is  silver-plated 
invar  so  the  relevant  work  function  is  that  of  silver,  which  is  4.26eV. 

Both  of  Spectrolab’s  arrays,  #8  and  #11,  will  be  biased.  The  first  array  is  their  WTC 
design,  which  is  expected  to  operate  better  than  conventional  cell  designs  at  high  voltages. 
The  4cmx4cm  cells,  which  are  177/xm  (7mil)  thick,  are  covered  by  a  152/im  (6mil)  thick 
coverglass  and  mounted  on  a  Ge  substrate.  The  second  array  contains  a  similar  cell  of 
conventional  design.  The  dimensions  of  both  the  cell  and  coverglass  are  the  same,  and  the 
adhesive  is  estimated  to  be  76//m  (3mil)  thick.  Like  the  ASEC  GaAs/Ge  cells,  these  cells 
also  have  a  silver-plated  invar  interconnector. 

The  Advanced  Photovoltaic  Solar  Array  (APS A)  is  the  final  high  voltage  biased  array 
studied.  These  thin-film  silicon  cells  are  expected  to  be  the  cells  of  the  future,  based  on  the 
mass  cost  savings  they  yield.  Built  by  TRW,  they  are  2.6cm  x  5. 1cm  of  only  56/im  (2.2mil) 
thickness.  More  importantly  for  arcing  considerations,  the  cell  CMX  coverglass  thickness 
is  merely  51^m  (2mil)  thick  and  the  DC  93500  cover  adhesive  is  estimated  to  be  76^m 
(3mii)  thick.  The  interconnector  is  again  silver-plated  invar. 

The  voltage  biasing  will  be  performed  in  pre-set  sequences  from  -500V  to  -I-500V,  with 
each  bias  step  held  for  23  seconds.  In  this  work,  predictions  were  made  at  25  volt  increments 
with  a  twenty  second  experiment  time. 

4.3  Predictions 

To  predict  the  arcing  rate  of  the  biased  arrays,  the  analytic  arc  rate  method  described  in 
Section  2.1.2  is  used.  Table  2  shows  the  given  and  assumed  cell  properties  necessary  for 
the  arc  rate  calculations.  To  determine  Cdieu  and  the  offset  of  due  to  the  surface 
charge  density  (see  Section  3.4),  the  numerical  model  must  be  run  initially.  Only  the  first 
scheme,  the  capacitance  matrix  method,  is  necessary  to  determine  the  polynomial  fit  for 
Cdxeu  ■  Since  the  surface  charge  density  is  not  incorporated  into  the  analytic  model,  an  offset 
of  Tcfecinumerical)  versus  Tcjcf {(analytical)  must  be  accounted  for.  The  case  chosen  by 
Cho^  closely  matches  the  numerical  and  analytical  Tf/et  values.  Since  Tf/pf  versus  dV 
curves  have  the  same  shape,  the  best  technique  is  to  determine  the  offset  of  each  curve 
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Table  2:  PASP  Plus  data  used  for  arc  rate  predictions 


Instrument  No. 

1,2 

3 

4,6 

8 

11 

36 

Cell  Type 

Si 

Si  WTC 

GaAs/Ge 

GaAs/Ge  WTC 

GaAs/Ge 

APSA 

Manufacturer 

RCA 

NASA 

ASEC 

Spectrolab 

Spectrolab 

TRW 

Cell  Size  (cm^) 

2x4 

8x8 

4x4 

4x4 

4x4 

2.6x5. 1 

No.  of  Cells 

20,  60 

4 

20,  12 

4 

8 

12 

Cell  Gap  i^m) 

500 

1000 

500 

1000 

500 

635 

^cell  (/^^) 

203 

203 

89 

178 

178 

56 

di  (Atm) 

153 

203 

102 

152 

152 

51 

d2  (Atm) 

37 

N/A 

51 

N/A 

76 

76 

3.5 

4 

4 

4 

4 

4 

2.7 

3 

2.7 

3 

2.7 

2.7 

')inax\ 

3.46 

4 

4 

4 

4 

4 

')max2 

3 

2 

3 

2 

3 

3 

Ernax,  (eV) 

330 

400 

400 

400 

400 

400 

Emax2  (fi^) 

300 

200 

300 

200 

300 

300 

(i>w  (eV) 

4.76 

4.85 

4.26 

4.5 

4.26 

4.26 

from  the  case  that  can  be  accurately  predicted  analytically.  Therefore,  the  entire  numerical 
model  is  used  at  a  specified  voltage  with  the  orbit  integration  scheme  run  for  several  3 
values,  typically  five,  to  determine  this  offset. 

The  orbital  data  for  the  analytical  model  requires  a  different  approach.  Since  the  orbit  is 
neither  circular  nor  of  constant  altitude,  the  orbital  parameters  must  be  determined  for  each 
point  along  the  orbit  at  which  the  arc  rate  is  to  be  calculated.  The  International  Reference 
Ionosphere  (DU)  program,  developed  at  Goddard  Space  Flight  Center  in  June  1987,  can 
calculate  the  ambient  density  n^,  the  electron  temperature  Tc,  and  the  ion  mass  rn^  if  given 
the  altitude,  month,  hour,  sunspot  number,  and  the  geodetic  coordinates.  Since  the  exact 
time  of  the  biasing  is  unknown,  the  parameters  are  averaged  over  the  entire  year  by  hour 
and  by  month.  The  times  that  gave  the  lowest  and  highest  ambient  densities  are  also  used 
to  define  the  lowest  and  highest  arcing  rates  which  might  be  observed.  The  variation  due  to 
different  sunspot  numbers  is  within  the  variation  due  to  ambient  densities.  To  determine  the 
coordinates,  as  weU  as  the  altitude  and  orbital  velocity,  a  separate  calculation  is  necessary. 
Orbital  mechanics  are  used  to  find  the  state  vector  in  inertial  coordinates  of  a  given  altitude. 
These  are  then  converted  to  the  earth-fixed  coordinates.  By  using  the  orbit  inclination. 
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longitude  of  the  ascending  node  and  the  orbital  eccentricity,  the  equations  of  motion  can 
be  integrated  using  a  fourth-order  Runge-Kutta  method  to  obtain  other  positions  along  the 
orbit. 

At  altitudes  higher  than  1000km,  the  IRI  model  is  useless  for  obtaining  particle  densities. 
At  these  altitudes,  however,  the  following  formula  developed  by  Al’Pert’  has  been  shown 
to  agree  well  with  experimental  measurements: 


Ueiz)  =  exp 


2kTr 


1  -h 


Ea 

z 


1  -h 


-f  In 


1 


z 


1  -h 


(5:t) 


In  this  equation  Ro  is  the  Earth’s  mean  radius,  is  the  electron  density  at  1000km,  go 
is  the  gravitational  acceleration  at  1000km,  Mo  is  the  mass  of  hydrogen  at  1000km,  is 
the  electron  temperature  at  1000km,  and  is  the  altitude  of  1000km.  Finally,  the  arc  rate 
program  requires  the  input  of  a  range  of  electron  densities.  Since  each  voltage  bias  will  be 
held  for  23sec,  the  electron  densities  at  either  end  of  that  range  along  the  orbit  are  used  to 
define  the  maximum  and  minimum  densities  expected. 

For  each  array  arc  rates  are  predicted  for  voltages  up  to  -500V  in  25V  increments  at 
intervals  of  600sec  around  the  orbit.  Near  apogee,  the  spacecraft  slows  and  the  arcing 
rales  have  less  variation  due  to  the  lower  ambient  densities  so  the  intervals  are  increased  to 
1200sec. 

The  arc  rate  results  for  the  conventional  arrays  are  presented  in  Figures  48-59.  For  each 
array,  the  first  figure  shows  representative  arc  rates  around  the  spacecraft  orbit  for  the  entire 
voltage  range.  Standard  deviation  error  bars  included  on  these  figures  are  based  on  typically 
100  calculations  for  each  voltage.  The  arc  rates  at  lower  voltage  biases  are  identical  since 
Tc fee  is  dominant.  The  diverging  curves  at  higher  voltages  indicate  non  is  beginning  to 
affect  the  arc  rates.  To  distinguish  the  arc  rates  at  these  higher  voltages,  a  second  figure  for 
each  array  shows  all  of  the  arc  rate  curves  in  the  range.  As  shown  in  Figure  48,  the  arc- 
rates  for  the  silicon  conventional  array  with  20  cells  range  from  an  average  of  0.004sec~  ‘  at 
200V  to  21  lsec“'  at  500V.  At  500V  bias,  the  arcing  rates  are  as  low  as  8.75sec'',  which  is 
calculated  at  apogee.  The  addition  of  40  more  solar  cells  to  the  array  increases  the  number 
of  emission  sites  where  arcing  can  occur,  thus  decreasing  the  time  between  arcs.  This  is 
clearly  seen  in  Figure  50,  which  shows  that  arcing  begins  at  voltages  as  low  as  175V  for 
array  #2.  The  highest  arcing  rate  calculated  is  366sec~'  at  perigee  for  a  voltage  bias  of 
500V,  while  at  apogee  the  arcing  rate  is  33sec“’.  For  ASEC’s  GaAs/Ge  array  of  20  cells 
(#4),  arcing  is  predicted  at  125V,  as  shown  in  Figure  52.  Accordingly,  the  arcing  rates  are 
much  higher  at  500V,  with  a  maximum  of  1382sec"'  and  a  minimum  of  50sec~^  The 
slightly  smaller  array  of  only  12  cells  (#6)  is  also  predicted  to  experience  arcing  at  125V, 
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but  its  arcing  rates  at  5(X)V  range  from  918sec-'  to  only  33sec"'.  The  arc  rate  predictions, 
shown  in  Figure  56,  for  Spectrolab’s  GaAs/Ge  array  of  8  cells  are  significantly  lower  than 
predictions  for  ASEC’s  GaAs/Ge  array.  Arcing  is  not  predicted  until  175V,  and  arc  rates  at 
500V  are  at  most  69sec“'  (at  perigee)  and  at  least  15sec-^  (at  apogee).  The  thin  APS  A  is 
predicted  to  experience  arcing  at  the  low  voltage  of  lOOV  in  the  lower  altitudes.  The  high 
arcing  rates,  however,  are  not  nearly  as  high  as  for  arrays  #4  and  #6.  The  lower  arc  rates 
are  much  higher,  indicating  a  higher  Te/ee,  but  the  arc  rates  saturate  around  300V  by  T^o„. 
Consequently  at  500V,  the  highest  arc  rate  is  predicted  to  be  only  413sec'''  and  the  lowest 
is  predicted  to  be  15sec“'. 

To  compare  these  results,  Figure  60  shows  the  arc  rate  predictions  for  each  conventional 
array  at  perigee.  The  GaAs/Ge  arrays  and  APSA  only  differ  in  Tf/ef.  properties  in  dielectric 
thickness  (See  Table  2).  As  expected,  the  thin  APSA  has  the  highest  arc  rate  predictions  in  the 
Tf  /ff-dominanted  range,  with  only  127/xm  dielectric  thickness.  Spectrolab’s  GaAs/Ge  cells 
with  228/u.m  thick  dielectrics  have  the  lowest  predicted  arc  rates,  and  ASEC’s  GaAs/Ge  cells 
with  1 53/rm  thick  dielectrics  are  predicted  to  experience  arcing  of  intermediate  intensity.  All 
of  these  arrays  have  silver-plated  interconnectors,  in  contrast  with  the  Kovar  interconnectors 
of  the  silicon  cells.  From  the  arc  rate  sensitivity  study  of  different  interconnector  work 
functions  in  Chapter  3,  however,  the  difference  in  arc  rates  between  these  two  metals  is 
negligible  (see  Figure  42).  The  silicon  cells  also  have  190^m  thick  dielectrics,  which  is 
between  the  thicknesses  of  Spectrolab’s  GaAs/Ge  cells  and  ASEC’s  GaAs/Ge  cells.  The 
predicted  arc  rates  of  the  silicon  arrays  also  fall  between  the  arcing  rates  of  those  arrays,  so 
the  dielectric  thicknesses  must  be  the  discriminating  factor  among  all  of  the  arrays  for  the 
range  dominated  by  Te/ee-  For  the  range  dominated  by  the  major  difference  between 
the  arrays  is  the  cell  frontal  area,  Acci/.  According  to  Eq.  (6)  the  smaller  APSA  cells  increase 
Tiotj,  thus  decreasing  the  arc  rate  saturation  level  from  that  of  the  larger  cells. 

To  predict  the  arc  rates  of  the  WTC  cell  arrays,  an  analytical  model  similar  to  the  one 
for  the  conventional  cell  arrays  should  be  used.  This  model  has  not  yet  been  determined  so 
the  arc  rate  predictions  for  this  cell  design  are  left  for  future  work. 
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Figure  48:  Selected  arc  rate  predictions  with  standard  deviation  errors  for  Si  conventional  array  #1 


Figure  49:  Complete  arc  rate  predictions  in  the  differentiating  voltage  range  for  Si  array  #1 
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Figure  50:  Selected  arc  rate  predictions  with  standard  deviation  errors  for  Si  conventional  array  #2 
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Figure  5 1 :  Complete  arc  rate  predictions  in  the  differentiating  voltage  range  for  Si  array  #2 
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Figure  52:  Selected  arc  rate  predictions  with  standard  deviation  errors  for  GaAs/Ge  conventional  array  #4 


Figure  53:  Complete  arc  rate  predictions  in  the  differentiating  voltage  range  for  GaAs/Ge  array  #4 
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Figure  55:  Complete  arc  rate  predictions  in  the  differentiating  voltage 
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Figure  57:  Complete  arc  rate  predictions  in  the  differentiating  voltage  range  for  GaAs/Ge  array  #1 1 
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Figure  60:  Arc  rate  prediction  comparison  for  all  PASP  Plus  conventional  arrays  at  350km 
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Chapter  5 
Conclusions 


The  emergence  of  high  voltage  solar  arrays  has  been  driven  by  recent  space  system  designs 
which  require  high  power  levels.  Consequently,  only  recently  has  the  negative  interaction 
of  arcing  on  high  voltage  solar  arrays  in  low  earth  orbit  been  a  concern.  Since  arcing 
can  eventually  destroy  solar  cells,  it  is  imperative  to  investigate  methods  to  reduce  or  even 
eliminate  it  as  well  as  to  develop  methods  to  predict  it  for  existing  solar  cells  which  might 
be  used  on  future  high  voltage  solar  arrays.  This  research  has  addressed  both  of  these  needs. 

All  of  the  research  presented  in  this  report  has  been  based  on  the  numerical  and  analyt¬ 
ical  models  discussed  in  Chapter  2.  These  models  accounted  for  the  ambient  ion  charging 
and  enhanced  electron  field  emission  (EFEE)  charging  processes  which  were  previously 
determined  to  be  the  two  main  factors  leading  to  arc  initiation.  The  numerical  model  simu¬ 
lates  this  arc  initiation  process  by  integrating  the  particle  trajectories  self -consistently  with 
the  dielectric  surface  charge.  Electrons  emitted  from  the  conductor  that  hit  the  dielectric 
surface  released  secondary  electrons,  creating  a  positive  surface  charge  which  further  en¬ 
hances  the  electric  field  at  the  triple  junction.  With  an  EFEE  charging  site  on  the  conductor 
near  the  triple  junction,  the  exponential  dependence  of  the  EFEE  current  on  the  electric 
field  causes  the  electric  field  to  run  away.  This  gives  rise  to  arc  initiation.  The  analytical 
model  represents  the  time  to  arc  initiation  by  the  sum  of  the  ambient  ion  charging  time  and 
the  EFEE  charging  time.  The  numerical  model  was  also  modified  to  represent  the  wrap- 
through-contacl  geometry  using  the  same  charging  processes.  The  electric  field  runaway, 
signifying  arc  initiation,  occurs  at  the  triple  junction  including  the  semiconductor  and  the 
coverglass.  Two  different  dielectric  charging  behaviors  are  exhibited  for  relatively  low  and 
high  field  enhancement  factors. 

The  models  for  conventional  solar  cells  were  used  to  determine  methods  of  arcing  mit¬ 
igation  in  Chapter  3.  Based  on  the  analytical  model,  five  cell  properties  were  varied  to 
determine  their  effect  on  the  arc  rate:  the  coverglass  and  adhesive  dielectric  thickness,  the 
interconnector  work  function,  the  secondary  electron  yield,  the  ratio  of  the  dielectric  con- 
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slants,  and  the  length  of  the  coverglass  overhang  above  the  interconnector.  The  dielectric 
thickness  was  found  to  affect  the  arc  rate,  but  significant  reduction  could  only  be  attained 
by  large  thicknesses  which  would  be  costly  in  mass  addition.  The  interconnector  work 
function  did  not  affect  the  arc  rate  much  for  typical  metals,  though  very  high  work  functions 
decreased  the  arcing  rates  and  very  low  work  functions  increased  the  arcing  rates  at  lower 
voltages.  Secondary  electron  yields  below  unity  eliminated  arcing  since  the  electric  field 
could  not  run  away;  yields  just  above  unity,  however,  did  not  reduce  the  arc  rate  significantly. 
Coverglass  to  adhesive  dielectric  constant  ratios  of  less  than  unity  caused  the  surface  charge 
near  the  triple  junction  to  remain  negative  for  a  considerably  long  time,  thus  significantly 
reducing  the  arcing  rates.  Overhanging  the  coverglass  above  the  interconnector  also  caused 
the  dielectric  surface  to  charge  negatively  if  the  overhang  was  longer  than  the  critical  over¬ 
hang.  Overhangs  less  than  the  critical  overhang  only  acted  as  thicker  dielectrics  with  no 
negative  surface  charging.  Shorter  overhangs  did  not  affect  the  arc  rate  much,  but  longer 
overhangs  significantly  reduced  arcing  rates.  Therefore,  arcing  rates  can  be  best  reduced 
by  increasing  the  dielectric  thicknesses,  selecting  dielectric  materials  with  coverglass  to 
adhesive  dielectric  constant  ratios  near  unity,  and  overhanging  the  coverglass  above  the 
interconnector.  Dielectric  secondary  electron  yields  lower  than  unity  will  eliminate  arcing 
at  all  voltages. 

Arcing  rates  for  actual  solar  cells  were  predicted  in  Chapter  4  for  a  near  future  space 
experiment,  the  Photovoltaic  Array  Space  Power  Plus  Diagnostics  (PASP  Plus)  experiment. 
Using  the  numerical  and  analytical  models,  the  high-voltage-biased  conventional  cell  arrays 
were  studied.  The  most  significant  factors,  which  discriminated  the  arcing  rates  in  the  EFEE 
charging-dominated  voltage  range,  were  the  cell  thickness  and  number  of  cells  in  the  array. 
As  expected  from  the  sensitivity  scans,  the  cells  with  thinner  dielectrics  had  higher  arcing 
rates  at  the  low  voltages.  In  addition,  the  increased  number  of  EFEE  emission  sites  due  to 
a  higher  number  of  cells  increased  the  arcing  rates. 
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